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0. Introduction. 

In the present note we address the issue of singular Lagrangians in analytical mechanics. Deriving 
^ ■ Hamiltonian formulations of physical systems with singular Lagrangians was attempted by Dirac and 
On ' Bergmann [1] . The aim was to obtain Hamiltonian formulations of relativistic field theories although 

(N : 

Dirac formulated his theory in terms of finite dimensional geometry. Applying Dirac procedures to 
relativistic mechanical systems we find that in most cases the resulting Hamiltonian description con- 

■ tains less information than was available in the Lagrangian formulation. We propose a version of the 
^\ \ Legendre transformation without this defect. 

On ■ In a recent paper Cendra, Holm, Hoyle, and Marsden [2] express the opinion that Lagrangian 

systems and Hamiltonian systems offer different representations of the same object. The Legendre 
0^' transformation is the passage from one of these representations to the other. We agree with these 
concepts. We also agree with the statement that "one should do the Legendre transformation slowly 
and carefully when there are degeneracies" . We think that our Legendre transformation is slow and 
careful enough to provide the correct Hamiltonian representation of relativistic mechanical systems. 

We provide an almost complete although somewhat superficial review of the geometric background 
for analytical mechanics. Complete coordinate characterizations of all structures are provided. Intrinsic 
constructions of most of the object are given. A more rigorous version of this material is in preparation. 
Related material can be found in [10] [12] [13]. 

■ Much of the material was developed in collaboration with G. Marmo at Istituto Nazionale di Fisica 
Nucleare, Sezione di Napoli. 

1. Geometry of tangent and cotangent bundles. 

Let Q be a differential manifold of dimension m. We use a coordinate system or a chart 

■.X^{q^){x) = {q\x),... ,q^{x)). (1) 
Each individual coordinate is a function 

q^:Q^R. (2) 

We ignore the fact that the domain of a chart could be an open submanifold of Q and not all of Q. 
Let F be a differentiable function on Q. The function 

i^o (9«)-i:]R" ^ R (3) 
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is the coordinate expression of the function F. It is a function of the coordinates {q'^{q)) S of a 
point q G Q. We define partial derivatives 

*-^^"('') 

These partial derivatives are functions on Q. 

The tangent bundle of a manifold Q is a manifold TQ. There is a mapping 

tq:TQ^Q (5) 

called the tangent fibration. Tangent vectors (elements of TQ) are equivalence classes of curves in Q. 
Two curves -f.R Q and Y:R Q arc equivalent if 7'(0) = 7(0) and D(/ o 7')(0) = D(/ o 7)(0) 
for each function f:Q ^ R. The equivalence class of a curve j:R ^ Q will be denoted by t7(0). 
Coordinates 

■.v^{q\v),... ,q"^{v)Jq\v),... ,5q"^{v)) (6) 

are induced by coordinates (g") in Q. If 7 is a representative of a vector v, then q'^{v) = 5'*(7(0)) and 
5q^{v) = D(g'^ o 7)(0). The tangent fibration is defined by 

rQ(t7(0)) = 7(0). (7) 
Fibres of this fibration are vector spaces. We have operations 

+:JQ X JQ^JQ (8) 

and 

•:MxTQ^TQ (9) 

with coordinate representations 

{q\5q^){v^+V2) = {q"" {vi) , Sq\vi) + Sq\v2)) (10) 

and 

{q^,Sq^){k-v) = iq^{v),kSq\v)). (11) 
We denote by TQ x TQ the set 

{{vi,V2) e Jg X TQ; tq{vi) = tq{v2)} (12) 

Since representatives of vectors (curves in Q) can not be added the construction of linear operations 
in fibres of tq is somewhat indirect. Let v = t7(0), vi = t7i(0), and V2 = t72(0) be elements of the 
same fibre TqQ = TQ^{q). We write 

V = V1+V2 (13) 

if 

D(/ o 7)(0) = D(/ o 7i)(0) + D(/ o 72)(0) (14) 

for each function / on Q. We have defined a relation between three elements of a fibre TqQ. This 
relation will turn into a binary operation if we show that for each pair {vi,V2) € TqQ x TqQ there is 
an unique vector v e TqQ such that v = Vi + V2- The coordinate construction 

(g« o 7)(s) = {q-{vi) + {Sq^{vi) + 6q-{v2))s) (15) 
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of a representative 7 of w proves existence. Let v = t7(0) and v' = 17' (0) be in relations v = vi +V2 
and v' = Vi + V2 with vi = t7i(0) and V2 = t72(0). Then 

D(/ o 7')(0) = D(/ o 7)(0) = D(/ o 7i)(0) + D(/ o 72)(0) (16) 

for each function / on Q. It follows that 7' and 7 represent the same vector v' = v. This proves 
uniqueness. Let v = t7(0) and u = tA(0) be elements of TqQ and let fc be a number. We write 

V = ku (17) 

if 

D(/o7)(0) = A;D(/oA)(0) (18) 
for each function f on Q. The coordinate construction 

{q^o^){s) = iq^iu) + k6q^iu)s) (19) 

shows that for each fc € M and u € TqQ there is a vector v G TqQ such that v = ku. U v = t7(0) and 
v' = 17' (0) are two such vectors, then 

D(/ o 70(0) = D(/ o 7)(0) = fcD(/ o A)(0). (20) 

It follows that the vector v is unique. 

Each curve 7: M ^ Q has a tangent prolongation 

:s^t7(- + s)(0). (21) 

The curve 7(- + s) is the mapping 

7(- + s):K^Q 

:s'^7(s' + s) (22) 

The vector t7(s) is the vector tangent to 7 at 7(s). The coordinate description of the prolongation is 
given by 

{q^,6q^) otj ={q^o J, D{q^o^)). (23) 

A mapping X:Q^ TQ such that tq o X:Q ^ Q is the identity mapping is called a section of the 
fibration tq . A section of the tangent fibration is called a vector field. 
Let P be a differential manifold with coordinates 

(py.P-^W (24) 

For each differentiable mapping 

a:Q^P (25) 

we have the tangent mapping 

Ta: TQ TP. (26) 

If 7: R — > Q is a representative of a vector v G TQ, then a o 7: M — > P is a representative of the vector 
Ta{v) e TP: 

Ta(t7(0)) =t(ao7)(0). (27) 
The coordinate definition of the tangent mapping is given by 

{p\ 6p^) o Ta = {a' o tq, {d^a^ o TQ)5q'') (28) 
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with a' = o a or by a simplified formula 



{p\ V) o Ta = {a\d^a^q''). (29) 
Einstein's summation convention is used. The commutative diagram 

(30) 




is a vector fibration morphism. 

A difFerentiable mapping a-.T ^ Q is called an immersion if at each point t G T the linear mapping 
Tjcj: TtT — » T(j(j)(5 obtained by restricting the mapping Ta to the fibre T^T = T^^(t) is injective. If 

(f):T^M'= (31) 

are coordinates in T and a'^ = q*^ o a, then a is an immersion if the matrix {dicr'^) is of maximal rank 
k. The image S = im(cr) C Q is called an (immersed) submanifold of Q of dimension k. A submanifold 
5* C Q is frequently given as a set 

S={q€Q;y^FA{q)=0}, (32) 

where Fa are to — A: functions on Q such that the matrix {di^FA) is of maximal rank m—k at points of S. 
A set 5 specified in this way is called an embedded submanifold. Submanifolds are usually assumed to 
be embedded. We will adopt the standard practice of not distinguishing elements of geometric spaces 
from their coordinates. Functions defined on these geometric spaces will be considered functions of 
coordinates. Instead of writing a formula (32) we will say that S satisfies equations FA{q'^) = 0. The 
tangent set of a subset S C Q (not necessarily a submanifold) is a subset of TQ. A vector v is in 
T5 if there is a curve 7: R ^ Q such that v = t7(0) and 7(5) € S for each s in a neighbourhood of 
€ M. We have tq{TS) = S. US is the image of an immersion a:T ^ Q, then TS is the image of 
Ta: TT — > TQ. The coordinates {q'^, Sq^) of elements of TS are related to coordinates {t^, SP) by 

q"" = a'^if), 5q^ = dja^{e)5t^ . (33) 

If S satisfies equations FA{q'^) = 0, then TS* satisfies equations d^FASq'^ = in addition to FA{q'^) = 0. 
A 0-form on Q is a function on Q. A 1-form on Q is a mapping 

A-.TQ-^R 

■.v^ {A, v) (34) 
linear on fibres of tq . The product of a 0-form F with a 1-form ^4 is a 1-form FA defined by 

{FA,v) = FiTQ{v)){A,v). (35) 
The differential dF of a function on Q is 1-form defined by 

(dF,t7(0)) =D(i^o7)(0). (36) 

The differential of the product FG of two functions is the 1-form FdG + GdF. Coordinates (Sq'^) 
in TQ are 1-forms. They are the differentials {dq'^) of coordinates (g") in Q. Each 1-form A can be 
expressed as a combination 

A = A^dq'' (37) 
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of these differentials. The coefficients A^, are 0- forms obtained from 

{A,v)^A^{v)dq^{v) (38) 
for each v G TQ. The differential of a function F{q'^) is the 1-form 

dF = dxF{q'')dq^. (39) 



A 2-form on Q is a function 



which is antisymmetric: 



B:TQ X TQ M 
:{vi,V2) 1-^ {B,vi AV2), (40) 

{B,viAv2) + {B,V2Avi)=0 (41) 



and linear in its first argument: 

{B, {kvi + k'v[) A V2) = k{B,vi A V2) + k'{B,v[Av2). (42) 

Linearity in the first argument and antisymmetry imply linearity in the second argument. The product 
of 0-form with a 2-form is a 2-form. The exterior product of 1-forms A^ and A^ is a 2-form A^ A A^ 
defined by 

{A' AA'',viAv2) = {A\vi){A^V2) - {A\v2){A^ ,vi). (43) 



Each 2-form B is a combination 



B = ^B«Adg« A dq\ (44) 



The coefficients B^x are 0-forms characterized by 

{B, V, A V2) = lB,xiSq^{vi)Sq\v2) - dq^iv2)5q\v,)) (45) 

and 

B.x + B^, = 0. (46) 

The exterior differential of a 1-form A is a 2-form dA. In order to construct the exterior differential 
we associate with each pair (ui, f2) € TQ x TQ a mapping %: ^ Q such that vi = tx{-, 0) and 

V2 = tx(0, •). The coordinate construction 

X^isi, S2) = q''{x{si,S2)) = q^ivi) + Sq'^{vi)si + dq''{v2)s2 (47) 
proves the existence of such mappings. We define curves 

^i:M^TQ 

:s^tx(-,s)(0) (48) 

and 

6:M^TQ 

■.s^tx{s,-m (49) 

with coordinate representations 

{^ns2),S^^{s2)) = (g'^(a(s2)),<5g\a(s2))) = {x^iO,S2),d,,x\0,S2)) (50) 



and 

msi),SC^{s^)) = {q^{Usi)),Sq\a{si))) = ix^{suO),d,,x\suO)). (51) 
For the mapping defined in (47) we have 

i^iis),S^^is)) = {q'^M+Sq-iv2)s,dq\v,)) (52) 

and 

(^2 (s), S^^{s)) = {q-{v^) + 6q-{vi)s, Sq\v2)) (53) 
The exterior differential is defined by 

(d^, vi A V2) = D(A, ^2)(0) - B{A, 6)(0) (54) 

Relations 

A^AA^ + A^AA^= 0, (55) 

d{FA) = dFAA + FdA, (56) 

and 

ddF = (57) 

are easily established for an arbitrary O-form F and arbitrary 1-forms A, A^, and A"^. The exterior 
differential of a 1-form A = Axdq^ is the 2-form 

d^ = dAx A dg^ = d^Axdq^ A dq^ = ^{d^A^ - dxA,)dq'' A dq\ (58) 

A 2-form which is the differential of a 1-form is said to be exact. 

A 1-form A is said to be dosed if dA = 0. If A is closed, then there is a neighbourhood V of each 
point qo and a O-form F on V such that = dF. This is as a consequence of the Poincare lemma. 

Let P be a differential manifold with coordinates (p^) and let a: Q P be a differentiable mapping. 
Let a* = o a. The pull back of a O-form P on P is the O-form P o a on Q. The pull back of a 1-form 
A on P is the 1-form a*AonQ defined by 

{a*A,v) = {A,Ta{v)). (59) 

If 

A = Aidp\ (60) 

then 

a*A = Aid^a'dq" (61) 
The pull back of a 2-form B on P is the 2-form a*B on Q defined by 



(a*P, vi A V2) = {B, Ta(wi) A Ta{vi)}. (62) 



If 

then 

The relations 
and 



P=-Pyd/Ady, (63) 



a*B = ^B.jd^a'dxa^dq'' A dq^. (64) 

d(a*P) = a*dP (65) 

d{a*A) = a*dA (66) 
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hold for a 0-form F and a 1-form A. Let C C Q be a submanifold. The mapping 

■.q^ q (67) 

is the canonical injection. The pull backs t^-^' '•c^' ^"^"^ ^C-^ ^-"^^ denoted by F\C, A|C, and i?|C 
respectively. 

The cotangent bundle of a manifold Q is a manifold T^Q. The cotangent fibration 

7Tq:T*Q^Q (68) 
is the vector fibration dual to the tangent fibration tq. The canonical pairing is a bilinear mapping 

(,):T*g X TQ^R 

■{f,v)^{f,v) (69) 

defined on the set 

T*Q X TO = {(/,t;)eT*OxTO; 7rQ(/) = rQ(^;)} (70) 

Each covector f e T*(5 = T^g^iq) is the differential dF(5) of a function — > R. Differentials 
(dg'*(g')) form a basis of the' vector space T*Q. Let (6^(9)) be the basis of the vector space TqQ dual 
to the base {dq'^{q)) in the sense that 

{dq-{q),ex{q))=S\. (71) 

Coordinates 

(3",/a):T*Q^R2'" ^72) 

are defined by 

('Z^/A)(/) - (9"(^Q(/)),(/,eA(7rQ(/)))) (73) 
The canonical pairing has the coordinate expression 

{f,v)=fx{f)Sq\v). (74) 

For the tangent bundle TT*Q of the cotangent bundle T*Q we have the tangent fibration 

r-r*Q:TT*g^T*Q (75) 

and the tangent mapping 

T7rQ:TT*Q^TQ (76) 
of the cotangent fibration ttq: T*(3 — > Q. The diagram 

Ttto 

TT*Q ^ Tg (77) 



is commutative. Hence, (rj*g(w), T7rQ(w)) G T*Q x TQ for each w G TT*(5. A canonical 1-form 
on T*(5, called the Liouville form, is defined by 

{'&Q,w) = {Tj^Q{w),TTrQ{w)). (78) 
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In the manifold TT*Q we have coordinates 

(g^PA,<5g^<5p,):TT*Q^M4n^ (79) 

related to the coordinates (9", /a) as the coordinates (q'^,Sq'^) in TQ are related to the coordinates 
(g") in Q. In terms o these coordinates, coordinates {q'^,fx) in T*Q, and coordinates {q'^,Sq^) in TQ 
we have the coordinate definitions of the fibrations Tj*q and Tttq: 

iQ^Jx)or^*Q = iQ^Jx) (80) 

and 

(g^V)oT7^Q = (5^V)• (81) 

It follows that 

{^Q,w) = Mw)6q^{w). (82) 

Hence, 

= f.dq\ (83) 

A 1-form ^ on Q is a function on TQ but it can be interpreted as a section A:Q — > T*Q of the 
cotangent fibration. In terms of this dual interpretation we state the following fundamental property 
of the Liouville form: 

A*^Q = A. (84) 

A manifold P and an exact, non degenerate 2-form uj form an (exact) symplectic manifold {P,u}). 
The 2-form w defines a mapping (3(p^^y.TP -^T*P characterized by the equality 

{fi{p,^)(u),v) = {uJ,u^v) (85) 

for vectors u G TP and v £ TP such that Tp{v) = Tp{u). The 2-form uj is said to be non degenerate if 
the mapping P(p^u>) is invertible. The cotangent bundle T*Q together with the 2-form 

ujQ = Mq = A dg« (86) 

form a symplectic manifold {T*Q,luq). In the cotangent bundle T*T*Q we use coordinates 

(g«,/A,a^,6^):T*T*g-^M4- (87^ 

induced by coordinates {q'^,f\) mT*Q. The coordinate definition of the mapping 

Vq,-«)^TT*Q^T*T*Q (88) 

is given by 

(g^/A,a^,6'')o/3(^*Q_^^) = (g«, /a, (5/^, -^g-^). (89) 
This mapping is invertible. Its inverse 

is defined by 

(g«, fx, Sq'^, 5f.) o /3(V*Q,.^) = A- (91) 

The Poisson bracket 

{F,G}:T*Q^M (92) 



of two functions F and G on T*Q is defined by 

{F,G}{f) = {dG{f),/3-^\^^^JdF{m. (93) 

It follows from the coordinate relation (91) that the Poisson bracket {F, G} of two functions F{q'^, fx) 
and G{q'^, fx) is the function 

dF^dG_ _ dG_dF_ 

or 

d^Fd'^G - d^Gd^F (95) 
with the symbol d'^ used to denote the partial derivative with respect to /„. 

2. Lagrangian submanifolds. 

A Lagrangian submanifold of a general symplectic manifold (P, lu) is a submanifold S* C P of di- 
mension dim(S') = i dim(P) such that u}\S — 0. This last condition means that the symplectic form 
oj evaluated on two vectors tangent to S vanishes. If S is the image of an immersion a:T ^ P, then 
U!\S = is equivalent to (t*uj = 0. 

A Lagrangian submanifold of {T*Q,ujq) is a submanifold S C T*Q of dimension m such that 
loq\S = 0. If S is the image of an immersion a:T ^ T*Q from a manifold T with coordinates {t") and 

(g^/A)oa=(a^aA), (96) 

then 

a*ujQ = d^a^d^a^dr A di^ = ^ {daa^df}a'^ - 5^a«a„a«) df A dt^. (97) 
If 5 is a Lagrangian submanifold, then the Lagrange brackets 

daa^df^a'' - dpa^do^a'' (98) 
vanish. Let f & S and let TfS C TfT*Q denote the space of vectors tangent to S at /. Let 

T}S = {a e T*T*g; \/^^^^^{a,w) = o} (99) 
be the polar of T/S". If m G TfS, then 

<Vq,-q)W'^) = {^Q,u^w)=0 (100) 



for each w e TfS. Hence, P(^t*q,ojq)0' f^) "'"/'^- Since dim(/3(j*Q ^^-)(T/5)) = dim(T/S') — m and 
dim(TJS) = dim(T*T*(5) — dim(T/S') = m, the spaces TfS and /3(x*q ujQ)0'f'^) equal. If F and 
G are functions on T*Q constant on S, then dF{f) and dG{f) are in Tj5 for each f G S. It follows 
that 

{P,G}|S' = 0. (101) 
If S is specified by equations Fa — 0, where Fa arc m independent functions on T*Q, then 

{F^,Fb}|5 = 0. (102) 

There are three categories of Lagrangian submanifolds of cotangent bundles generated by increasingly 
complex objects. 

I. Lagrangian submanifolds generated by functions. 
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Let f7 be a function on Q. The image S of the differential dU -.Q ^ T*Q isa, Lagrangian submanifold 
of {T*Q,ojq) since dini(S') = m and 

(dUfujQ = {dUfd'dQ = didUf-dq = ddU = 0. (103) 

The submanifold S is said to be generated by U. In terms of coordinates (g", f\) the set S is described 
by equations 

fx = dxU{q^), (104) 
equivalent to the simple version of the principle of virtual work 

fxdq'' = 5[/(g«) = aA[/(g«) V, (105) 

where the virtual displacements Sq^ are coordinates of a vector v £ TQ. 

Let S = im(o-) c T*Q be the image of a 1-form interpreted as a section a:Q ^ T^Q of the cotangent 
fi-bration. From 

{afojQ = {afd'dQ = d{(j)*i)Q = da (106) 

it follows that if S is a Lagrangian submanifold, then for each element fo & S there is a neighbourhood 
C T*(5 of /o and a function U on Q such that S (IW = im{dU) fl W. This is a version of the 

Poincarc lemma. 

II. Lagrangian submanifolds generated by constrained functions. 

Let C C Q be a submanifold of dimension k and let [/: C ^ M be a differentiable function. The set 

S={f€T*Q; g = 7rQ(/)eC,V,eT,ccT,Q(/'^) = (dC^'^)} W 

is an afBne subbundle of the cotangent bundle T*Q restricted to C. At each point q G C the fibre 
Sq = S m*Q is an affine subspace of T*Q modeled on the vector subspace T°C C T*Q oi dimension 
m — k. It follows that S' is a submanifold of T*Q of dimension to. We choose a function U:Q^M. 
such that U\C = U and define functions U = U ottq on T*Q and U = U\S on S. The fimction U does 
not depend on the choice of the function U, it can be defined directly by U{f) = U{'KQ{f)) for each 
f € S. If w e TS, then T-kq{w) e TC since ■kq{S) = C. From 

('(9q,w) = (T-pH<Q(!i7),T7rQ(w)) 
= {dU,lT,Q{w)) 
= {dU,lTZQ{w)) 

= {dU,w) 

= {dU,w) (108) 



it follows that 
and 



'&Q\S = dU (109) 



luqIS = ddqlS = d{-dQ\S) = ddU = 0. (110) 

Hence, 5 is a Lagrangian submanifold of {T*Q,coq). 

Given a function U{q'^) and m — k independent functions FA{q'^) such that the set C is described 
by the equations FA{q'^) = we write the principle of virtual work 

FA{q^) = 

fx6q^ = dxU{q^)5q^ (111) 
dxFA{q'')Sq^ = 
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for the set S. Coordinates {q'^,f\) of elements of S satisfy the variational principle with arbitrary 

virtual displacements Sq'^ satisfying the last equality. This last equality indicate that the virtual 
displacements are coordinates of vectors tangent to C. Using Lagrange multipliers A"^ we write the 
equations for S in the form 

= 

fx = dxU{q'') + dxFA{q'')X 



(112) 



Let (t*) be the coordinates in C and let q"^ = a'^{t'^) be the coordinate expression of the canonical 
injection of C in Q. If U{t) is the internal energy, then S is represented by 

q" = a^if) , , 

• 113 
MV(f) = a,[/(f). 

Let C C Q be a submanifold and let S be an afiine subbundle of the cotangent bundle T*Q restricted 

to C modeled on the vector subbundle T°C of T*(5 restricted to C . If S' is a Lagrangian submanifold 
of (T*(5,a;Q), then iJqIS' is closed. Let /o be an element of S and let C T*(5 be a neighbourhood of 
/o and U a function on 5 fl W such that 'dQ\S(^W = dU. We choose the neighbourhood W to have a 
connected intersection Wq = SqCiW with the fibre Sq = S r]T*Q for each q in V = ttq^S DW). The 
restriction of t^g to the fibre T*Q is the zero form since (il^QjIu) = if TTrQ{w) = 0. Consequently, 

dU\Wq = '&Q\Wq = (114) 

and the function U is constant on the connected set Wq. This permits the introduction of a function 
U onC such that U{f) = U{-KQ{f)) for each / e 5 n VK. The set 

{/ e T*Q; q = WQ{f) e C, V,eT,CcT,Q(/'^') = mv)} (115) 

intersected with W is the intersection of S with W. We have obtained an extension of the Poincare 
lemma to constrained Lagrangian submanifolds. 

III. Lagreingian subrriEinifolds generated by Morse families. 

Let r]:Y ^ Q he a, differential fibration with coordinates {q'^, y"^) adapted in the sense that 

{qn°r,= {qn. (116) 

where the coordinates {q'^) on the right hand side are coordinates in Y. Let U:Y ~> M.he a function 
interpreted as a family of functions defined on fibres of the fibration rj. The family is called a Morse 
family if the A; x (m + fc) matrix 

8^U O^U ^ ^^^^^ 



dy^dy^ dy^dq'^ 
is of maximal rank. A Morse family generates a set 

5 = {/ e T*Q; 3,,v,,„,V,,t„f(/'T»7(^)) = {dU,z)} . (118) 

The critical set 

Cv{U, ri) = [ye Y; V^eT„yTr;(u;) = ^ {dU, w) = o} (119) 
of the Morse family is a submanifold of Y of dimension m. A mapping 

K:CT{U,r])^T*Q (120) 
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such that 'KQ{K{y)) = r}{y) is defined by 

{K{y),v) = {dU,w), (121) 

where v is any vector in T^iy) and w G Tj,y such that Trfiw) = v. This mapping is an immersion and 
S = im(K). Let y G Ct{U, rj) and w GTy Cr(?7, rj). From 

= {Tj*Q(TK{w)),TTrQ{TK{w))) 
= {K{y),Tr]{w)) 

= {dU,w) (122) 

it follows that 

K*dQ = dU\CT{U,ri). (123) 
The set S is an immersed Lagrangian submanifold of {T*Q,lvq) since 

K*LJQ = (124) 

and dim(5) = dim{Ci {11,7])) = m. 

It follows from a theorem of Hormander [4] [7] that for each clement /o of a Lagrangian submanifold 
of {T*Q, u)q) there is a neighbourhood W C T*Q and a Morse family U:Y ^ R of functions on fibres 
of a fibration r]:Y ^ Q such that 5 and the Lagrangian submanifold generated by U coincide in W. 
This is an extension of the Poincare lemma. 

The coordinates {q'^, fx) of elements of S satisfy equations 

fx = dxU{q^,y^) 
= 5sC/(9^/) (125) 

derived from the variational principle of virtual work 

fxdq^ = 5U{q'\ y^) = dxU{q\ y^)Sq^ + dsUiq^, y^)5y^ (126) 

with some values of the variables (y^) and all variations (Sq^'^jSy^). The symbol Oa stands for the 
partial derivative with respect to y^. Equations (125) imply the equalities 

dfx = d^dxU{q\.y^)dq'' + dBdxU{q\ y^)dt/^ 

= a;,aBC/(g^J/^). (127) 

Consequently, 

OjQ\S = dfxhdq^\S = d^dxU{q^,y^)dq^^ hdq^ =0. (128) 

It follows from the maximality of the rank of the matrix (117) that dim(S') = to. 
Note that the afRne subbundle (107) is generated by the Morse family 

U{q^,y^) = U{q^) + FA{q^)y^. (129) 

The rank of the matrix 

^ = fo f^) (130) 



QyAQyB dy^dq'^ J \ dq'^ 

is maximal due to independence of the functions F^iq'^)- The function (129) depends linearly on 
the unrestricted variables (y^). This is the characteristic feature of a Morse family equivalent to a 
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constrained generating function. There is little difference between the variables (y^) and the Lagrange 
multipliers (A'^). 

A Morse family generating a Lagrangian submanifold is not unique. It is frequently possible to 
reduce the dimension of the fibration 77. Reductions are based on the following observation. Let S be 
a Lagrangian submanifold of {T*Q,ujq) generated by a Morse family U:Y ^ M. oi functions defined 
on fibres of a fibration 77: F — + Q. If the critical set Ci{U, rj) is the image of a section a:Q^Yoir], 
then S is generated by the function [/ o cr: Q — > M. If / e 5 and q = ngif), then / = K{a{q)) and 

(/, v) = («(a(g)), v) = {dU, Ta{v)) = {d{U o a), v) (131) 

for each v G TqQ. Hence, / = d{U o (7){q). This shows that S C im(L'' oa). If / = d(L'' o a){q), then 

(/, v) = (d(C/ o a), v) = (dC/, Taiv)) = {K{a{q)), v) (132) 

for each v (E TqQ. Hence, / = K{(T{q)). It follows that im{U oa) C S. It may happen that the fibration 
77 is the composition rj" o rj' of fibrations rj': Y ^ Y' and rj": Y' ^ Q and that the critical set Cr(C/, 77') 
is the image of a section cr:Y' ^ Y of t]' . In this case the Lagrangian submanifold S is generated by 
the Morse family {/ o cr: y — > M of functions on fibres of rj". 

3. Statics of mechanical systems. 

Let Q be the configuration space of a static mechanical system. Elements of the cotangent bundle 
T*Q are the generalized forces applied to the system. The constitutive set of a static system is subset S 
(usually a submanifold) of the cotangent bundle. An element / e S* is the generalized force which when 
applied by an external controlling device will maintain the system in equilibrium at the configuration 
q = 7rQ(/). The constitutive set provides a complete characterization of the response of the static 
system to external control represented by generalized forces applied to it. The knowledge of equilibrium 
configurations of an isolated system does not characterize the system completely. Two systems may 
have the same equilibrium configurations and yet respond differently to external control. 

The system is said to be reciprocal if i-uq\S = 0. Let wi and W2 be vectors tangent to S such that 
Tj*q{w2) = Tj^giwi). Let Siq'^ = 5q'''{wi), Sif^ = df^iwi), 52q'^ = 6q'^{w2), and 62/^ = <5/k(w2). 
The equality 

Sifj2q'' = 62fJiq^ (133) 

derived from {uiq,wi A w;2) = is a reciprocity relation. The system is said to be potential if S is 

a Lagrangian submanifold generated globally by a generating function, a constrained function or a 
Morse family. The generating function is interpreted as the internal energy of the system. A potential 
system is reciprocal. 

In the following three examples the configuration space is an affine Euclidean plane with Cartesian 
coordinates {x,y). Coordinates {x,y,f,g) are used in T*Q. 

Example 1. The function 

U{x,y) = ^{x^+y^) (134) 

is the internal energy of an elastically suspended material point. The constitutive set S is the La- 
grangian submanifold generated by U. It is described by equations 

f = kx, g = ky (135) 

derived from the principle of virtual work 

fSx + gSy = dxU{x, y)5x + dyU {x, y)5y. (136) 

A 
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Example 2. Let C c Q he the circle 

a;2+y2 = a^. (137) 

Let _ 

U{x,y) = ky (138) 

represent the internal energy of a material point constrained to the circle 

x'^+y^ = c?. (139) 



From the variational principle 



+ 2/^ = 

/fo + g^y = k5y (140) 
x5x + ySy = 



we derive equations 



x^+y^ = 

f = \x (141) 
g = k + Xy 

for the constitutive set S with a Lagrange multiplier A. With the parametric representation 

X = a cos 1? 

(142) 

y = a sin 1? 

we obtain the expression U {{)) — ka sin i9 for the internal energy and the variational principle 

X = a cos d 

y = asini} (143) 
— fa sin '&5'& + ga cos '&S'& = ka cos iJ^t? 

equivalent to (140). The constitutive set is generated by the Morse family 

U{x,y,X) = ky + ^{x^+y^-a^). (144) 

▲ 

Example 3. The function 

U{x,y,^) = ^{{x - acos^f + {y - asin^f) (145) 

is the internal energy of a material point tied elastically to a point left to move freely on the circle 

a; = acosi?, j/ = asini?. (146) 

The function f7 is a Morse family of functions of the variable since the rank of the 1x3 matrix 

/ d'^U d'^U d'^U \ 

'Mdx j = (Ma;cost? + j/sint?), kasind, -kacosd) (147) 

is 1. From the principle of virtual work 

f5x + gSy = 6U{x,y, •&) = k{x — acos'&)5x + k{y — a) sm'&)dy + ka{x sm-d — y cosi?)5i? (148) 
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we obtain equations 

/ = k{x — acos'd) 

g = k{y - asvcid) (149) 
= kaixsin^ — ycos^) 

for the constitutive set S. Equations 

X = p cos 'd 
y = psm.'d 

f If ^ , ^^50) 
/ = k[p — a) cos w 

g = k{p — a) sint? 

represent a mapping a from M? to T*Q. The set S is the image of this mapping. The matrix 



dx dy df dg 

'dp dp 'dp 'dp \ = ( «in79 kcos-& fcsini? 

dx dy df dg j psini? pcos-d —k{p — a)sm'& k{p — a) con d 

m m d^ M 



(151) 



is of rank 2. This indicates that S is an immersed submanifold. With the exclusion of points corre- 
sponding to p = the set S is the union of images of two sections of ttq corresponding to the two 
different signs in the formulae 

(152) 

\/x^ + y2 V / 
With the exclusion of points corresponding to x^ + y'^ > a^, S is the set of points satisfying equations 



f2 . 

F^ix, y,f,g) = y- f {VTT? - ka) = 0. 



FS{x,y,f,g) = x- J _ (x//^ + - fca) = 

''V/ +f (153) 



The functions F° and F° are obviously independent. It follows that S is an embedded submanifold. 
The rank of the Jacobian 

/ dx dy 

'd'p 'd'p \ _ ( cosi} smy , ^^^^^ 



dx dy_ j \^-psint? pcost? 
V dp dp 

of the mapping ttq o a represented by 

X = p cos "& 

. , (155) 
y = psmxf 

changes from 2 to 1 at p = 0. This indicates the presence of a Lagrangian singularity above the point 
with coordinates {x,y) = (0,0). A 



4. Differential equations. 

The tangent fibration Tq'.TQ Q was introduced in Section 1. Elements of the tangent bundle 
TQ were interpreted as virtual displacements. We return to the topic of tangent vectors this time 
interpreted as velocities. Coordinates (g**, 5q^): TQ — > M^™ introduced in Section 1 will be now denoted 
by {q'^,q^). 
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A first order differential equation in Q is a submanifold D c TQ. A curve j:I—^Q defined on an 
open interval / C M is said to be a solution of D if for each t E I the vector tjit) tangent to 7 at 7(f) is 
an element of D. A differential equation D is said to be integrable if for each v £ D there is a solution 
7: 7 — > Q of such that v = t7(to) for some to S /. 

Not all differential equations are integrable. Let D C TQ be a differential equation and let C be the 
set tq{D). If V G D and D is integrable, then there is a solution 7: / ^ Q of £> such that v = t7(fo) 
for some to G I. Since t'y{t) G D for each t G I, it follows that 7(f) e C for each t G I. Consequently 
t'j{t) E TC for each t E I and v = t7(io) E TC. We have shown that the condition D C TC is 
necessary for integrability of the equation D. This condition is sufficient for a class of differential 
equations described below. 

The image D = im{X) of a vector field X: Q — > TQ is an integrable differential equation. Let C C Q 
be a submanifold and let D be the union 

y {im(X,|C)} (156) 

of a family of vector fields 

Xa-.Q^TQ (157) 
restricted to C. If D c TC, then each field induces a vector field 

Xa-C^TC 

■.q^X^iq) (158) 
since im(X„|C) C D C TC. The differential equation im{Xa) is integrable for each a G A and 

D=\J {im(X„)}. (159) 

Hence, D is integrable. 

If the necessary condition D C Ttq{D) is not satisfied, then the reduced equation D D Ttq{D) is 
closer to being integrable although the condition D n Ttq{D) C Ttq{D n Ttq{D)) is not necessarily 
satisfied. This observation suggests the following algorithm for extracting the integrable part of a 
differential equation. We consider the sequence of sets 

C° = tq{D),c' = tq{D n TC°), ...^C" = tq{D n TC'^"'), . . . (160) 

and the sequence of differential equations 

= D,D^ = Dn TC", ...,d'' = Dn TC''~^, . . . (161) 

It may happen that after a finite number of steps the sets in the sequence (160) are all equal to a set 
C. This set satisfies the equality 

C = TQ{DnTC). (162) 

If the differential equation D = D H TC is integrable, then it is the integrable part of D. 

In Section 11 we give an example of a version of the above algorithm applied to a Hamiltonian 
system. 

Other algorithms for extracting the integrable part of a differential equations have been designed. 
They require the use of higher order tangent vectors. 

The second tangent bundle of a manifold Q is the set T^Q of equivalence classes of curves in Q. 
Two curves -f.R ^ Q and j'-.R^ Q are equivalent if 7'(0) = 7(0), D(/ o 7')(0) = D(/ o 7)(0), and 
D^(/ o 7')(0) = D^(/ o 7)(0) for each function f:Q^R. We use coordinates 
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in T^Q. If 7 is a representative of a second tangent vector a G T'^Q, then (/"(a) = q'^{'y{Oj), q^{a) = 
T){q'^ o 7)(0), and (j'^(a) = D^((j^ o 7)(0). The equivalence class of a curve 7: K — > Q will be denoted by 
t^7(0). Each curve 7: R ^ Q has a second tangent prolongation 

■.t^t^^{-+t){0) (164) 
The coordinate description of the prolongation is given by 

{q\ q\ r) o t"^ = {q'^ o 7, D(g^ o 7), D^q'' o 7)). (165) 
The second tangent fibration is the mapping 

:t7(0) 7(0). (166) 

There is also the fibration 

t^2q:T^Q^TQ 

:t^7(0) h^t7(0). (167) 

A second order differential equation in Q is a submanifold E C T^Q. A solution is a curve ^: I ^ Q 
such that t^j{t) e E for each t in the open interval / C M. The concept of integrability is easily 
extended to second order equations. The image im{X) of a section 

X:TQ^T'^Q (168) 

of the fibration t^2Q is an intcgrable differential equation. 

Elements of the iterated tangent bundle TTQ arc equivalence classes of curves in TQ. Coordinates 

(5'=,5^,g''',g"'):TQ-^M4m (^gg) 

will be used. These coordinates are related to coordinates {q'^,q^) as coordinates {q'^,q^) are related 
to coordinates {q'^). We have fibrations 

TTQ-.JTQ^TQ (170) 
TrQ:TTQ-*Tg (171) 

(9^g^)°TTQ = («^«') (172) 
(9^a^)oTTQ = (g^g'^). (173) 



and 

with coordinate representations 
and 



There is an useful immersion Ag of T^Q in TTQ. This immersion associates with a second tangent 
vector a = t^7(0) the vector w = tt7(0) tangent to the prolongation t7 of the curve 7 at t7(0). The 
formal definition is expressed in 

Aq:T2Q^TTQ 

:t27(0) tt7(0). (174) 

Prom rrQ(tt7(0)) = 17(0), TTQ(tt7(0)) = t7(0), and T^2Q{t^l{0)) = tj{0) it follows that 

n-QoAQ = TtqoAq =tSq. (175) 
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Let D c TQ be a differential equation. The set 

PD = Aq1(T£))cT20 (176) 

is a second order differential equation called the prolongation of D. If the differential equation is given 
in the form 

D = {v& TQ; y.fi{v) = 0} , (177) 
where are functions on TQ, then TD is the set 

{w e TTQ; yjiiTTQiw)) = 0,a^/i(TTQ(«;))g'^(a) + d^MTTQ{w))q"^{a) = 0} (178) 

and PD is the set 

{a e T^Q; yj.ir^Qiaj) = 0,d^MT\Q{a))r{a) + di^f^iT^Q^aWia) = 0} . (179) 

The symbol 5^ stands for the partial derivative with respect to g**. 
The inclusion 

tSq(P£>)c£> (180) 

follows from 

t\q{9D) = tjq{\q{PD)) = tjq{\q{Xq\1D))) c tjq{1D) = D. (181) 

Let D C TQ be an integrable equation. If u S D, then there is a solution 7: / ^ Q of Z? such 
that t7(0) = V. We have \Q{\?-i{Q)) = tt7(0) G ID since t7(i) e for each f e /. It follows that 
V = t7(0) = T^2Q{t^l{<^)) e t12q(P£'). Hence, 

D^t\q{PD) (182) 

if is integrable. We have established a necessary condition 

t\q{PD) = D (183) 

for integrability of a differential equation D C TQ. If this condition is not satisfied, then the integrable 
part of Z) is a subset of the set t^2q{PD). The set r^^qiPD) is a subset of D closer to the integrable 
part without being necessarily integrable. These observations suggest a new algorithm for extracting 
the integrable part of a differential equation. We introduce the sequence of differential equations 

& = D,D^= t\q{PD°), ...,&= t\q{P&-^), . . . (184) 

It may happen that after a finite number of steps the sets in the sequence (184) are all equal to a set 
D. It may happen that D is the integrable part of D. 

5. The iterated tangent bundle. 

We have already introduced the iterated tangent bundle TTQ and the coordinates 

{q'',q^,q'^,q'l-TQ^R'^'". (185) 

The fibration 

ttq:TTQ^TQ (186) 
is a vector fibration. We have the operations 

+:TTQ X TTQ ^ TTQ (187) 
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and 

• : K X TTQ ^ TTQ (188) 

with coordinate representations 

(9% q\ q'liwi + w^) = iq^{w,), q\w,), q'^{w{) + g'^^), n^i) + (189) 

and 



The diagram 



(9^ q\ q"', r){k ■ w) = {q^{w), q\w), kq"^{w), kr{w)). 



Ttq 

TTQ ^ TQ 



(190) 



(191) 



'n-Q 



TQ 



Q 



is a vector fibration morphism. 
We show that the mapping 

Ttq: TTQ ^TQ 
is a vector fibration by constructing operations 

+: TTQ X TTQ ^ TTQ 

(Ttq, Ttq) 

and 



(192) 
(193) 



•:MxTTQ^TTQ. (194) 

Let wi and W2 be elements of TTQ such that 'Jtq{w2) = 'Ttq{w\). It is possible to choose curves 
^i:R TQ and £,2'^ — > TQ such that Wi = t^i(O), W2 = t^2(0) and rg 0^2 = tq 0^1. The coordinate 
constructions 

(g", q^) o a = {q^iwi) + q'^{wi)s, q\wi) + q'\w,)s) (195) 

and 

(9", q^) o 6 = {q^{w2) + q'^{w2)s, q\w2) + q'\w2)s) (196) 
provide an example. The operation + is defined by 



u;i+«;2=t(a+6)(0). 



The operation ~ is defined by 



(197) 
(198) 



k~tm = mm- 

Coordinate representations of these operations are given by 

{q\ q\ g'^ + W2) = (q^iw^), q\wi) + q\w2), q'^{w^), riw,) + ^ (^^2)) (199) 

and 

iq^,q\q"',r){k~w) = {q^ (w) , kq\w) , q' ^ {w) , kr{w)) . (200) 

The diagram 



TTQ > TQ 



(201) 



TQ 



TQ 



TQ 



Q 
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is a vector fibration morphism. 

Elements of the iterated bundle TTQ are equivalence classes of curves in a set of equivalence classes 
of curves in Q. A simpler representation of these elements is needed. Let I^'^ ^ Q be a differentiable 
mapping. For each s e M we denote by t^^'^^xi^j^) the vector tx(s,-)(0) e TQ. For each t G M. we 
denote by t^^''^)x(0, i) the vector tx{-,t){0) G TQ. We have curves 



t(o>i) 



X(-,0):K^TQ 



(202) 



and 



t(i.o) 



X(0,-):M^TQ. 



(203) 

Vectors tt(°'i)x(-, 0)(0) e TTQ and tt(i'°)x(0, ^(O) e TTQ will be denoted by tt(°'i)x(0, 0) and 
tt^^'°)x(0, 0) respectively. For each w € TTQ there is a mapping x'-^"^ ^ Q such that w = tt^°'^^x(0) 0). 
The mapping specified by coordinate relations 

{q" o X)(s, t) = (q^iw) + r{w)t + q"'{w)s + q"'{w)st) (204) 

has the required property. We consider mappings X'^"^ ~^ Q and x'-^'^ ^ Q equivalent if 

tt(0'i)x'(0,0)=tt(°'i)x(0,0). (205) 



These mappings are equivalent if 



and 



x'(o,o) = x(0,0), 
D(i^o)x'(0,0) =d(i^o)x(0,0), 
D(0'1)x'(0,0)=d(o-i)x(0,0), 

D(i'i)x'(0,0) = D(i'i)x(0,0). 



(206) 
(207) 
(208) 

(209) 



We have obtained an efficient representation of elements of TTQ. In terms of this representation we 
define the canonical involution 



KQ-.TTQ^TTQ 

:tt(0'i)x(0,0) ^ tt(i'0)x(0,0) =tt(0'i)x(0,0), 



with 



X:R2_Q 
: (s,t) 1-^ x(^,s)- 

The coordinate expression of this involution is given by 

iq^,q\q'^r)oKQ = {q-,q'\q^^,r)- 

The commutative diagram 



TTQ 
TQ 



-^TTQ 

= TQ 



(210) 



(211) 



(212) 



(213) 



is a vector fibration isomorphism. The diagram 
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ITQ TIQ 

TQ TQ 

is the inverse isomorphism. For a differentiable mapping a:Q ^ P we have 

TTa(tt(°'i)x(0,0))=tt(°'i)(aox)(0,0) 

and 

Kp o TTa = TTq o kq. 
Let A be a 1-form on Q. A 0-form i^A on TQ is defined as the function 

iTA{v) = {A,v). 



(214) 



(215) 
(216) 

(217) 



Let -B be a 2-form on Q. If w € TTQ, then (rTQ(w),TrQ(w)) € TQ x TQ since tqoTtq = tqottq. 
A 1-form itB on TQ is defined by 



(iT-B,w) = {B,rTQ{w) ATtq{w)). 



(218) 



Let F, ^, and B = dA be a 0-form, a 1-form, and an exact 2-form on Q respectively. We define a 
0-form dxF, a 1-form d^A, and a 2-form dyS on TQ by 



drF = irdF, 
drA = irdA + dirA, 



and 



drB = dirB = dirdA = ddxA. 
The coordinate expression of the function dxF is 

dTF{q^,q^)=dxF{q^)q\ 

If 



and 



then 



and 



A = A^{q'^,p,)dq^ 



B=-B^x{q'',Pu)dq^Adq^, 



\tA = A^{q>',p^)q'', 
drA = dxA^q^dq", 
\TB = B^x{q'',PuWdq^, 



dirdA = -d^,{dKAx - d\A^)q^'dq'' A dq^ + {OkAx - dXA^)dq'' A dq^ 



Each 1-form on Q can be expressed as a sum of products FdG and from 

(dT(i^dG),tt(°>i)x(0,0)) = ^ (F{x{0,t))-^^G{x{s,t)) 



(219) 
(220) 

(221) 

(222) 

(223) 

(224) 

(225) 
(226) 
(227) 

(228) 
(229) 



|s=0,t=0 
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and 



it follows that 



Hence, 



for each l-form A. 



^^{FdG,txi;tm)\t=o = I (F(x(0,t))^G(x(s,t)) (230) 



(dT(FdG),tt("^i)x(0,0)) = ^(i^dG,tx(-,i)(0))|*^o- (231) 
(dTAtt(°'i)x(0,0)) = ±{A,tx{;tm)\t=o (232) 



6. A geometric framework for analytical mechanics. 

Let Q be a manifold of dimension m. We have already described the geometry of the tangent bundle 
TQ, the cotangent bundle T*Q and the tangent bundle JJ*Q of the cotangent bundle T*Q. The 
present section is devoted to the study of the canonical symplectic structure of the bundle TT*Q. We 
will use coordinates 

(g'^,g^):TQ^R2m^ (233) 
{q^,Pxy.T*Q^R'"', (234) 



and 



(<z^PA,^^P.):TrQ^M4™ (235) 



in the manifolds TQ, T Q and TT Q. The mappings Tj*q and Tttq have now the coordinate expres- 
sions 

(g^PA)or^-*Q = (g^PA) (236) 



and 

We introduce the exact 2-form 



(«^r)oT7^Q = (g^g^). (237) 

dxcoQ = diTWQ. (238) 

It will be shown that this 2-form is non degenerate. The manifold TT*Q with the form drwg form 
a symplectic manifold (TT*Q, d^t^g)- We believe that the symplectic form dxujQ is the only natural 
symplectic form in TT*Q. The discovery of a second symplectic structure in TT*Q was announced 
in a recent Springer- Verlag publication [9]. We have not been able to identify the second symplectic 
structure. We strongly suspect that this announcement is false. The formula 

de, = d^'' A dq" + dp^ A dp^ (239) 

for the Marsden-Ratiu symplectic form does not seem to have an intrinsic meaning since elementary 
rules of tensor calculus have been violated. We have the coordinate expressions 



and 



dr^Q = p^q^ + p^r (240) 



drujQ^dp^Adq'' +dp^Adq''. (241) 



The fibration Tj*g: TT Q T Q is a vector fibration. We will construct a vector fibration structure 
for the fibration Tttq: TT*Q — > TQ. For two vectors z\ G TT*Q and Z2 S TT^Q such that T'kq{z2) = 
TirQizi) it is possible to choose representatives Ci: K ^ T*(5 and C2: M ^ T*Q such that zi = tCi(O), 
Z2 = tC2(0) and TTQ o (2 = T^Q o Ci- An example is provided by the coordinate constructions 

(q^Px) o Ci = {q'^izi) + r{zi)s,px{zi) +Pa(^i)s) (242) 
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and 

An operation 

is defined by 
An operation 
is defined by 



{q^Px) 0C2 = {q^{z2)+r{z2)s,px{z2)+Mz2)s). (243) 



+:TT*Q X TT*0-*TT*0 (244) 

(TwqJvq) 



^i + ^2=t(Ci+C2)(0). (245) 
•:RxTTg^TTg. (246) 



A:-tC(0)=t(fcC)(0). (247) 
Coordinate representations of these operations are given by 

(?" , Pa , 5" , P..) (-21 + ^2 ) = (g" (-Z1 ) , Pa (-21 ) + PA {Z2 ) , 5" (-zi ) , P^ {zi)+Mz2)) (248) 

and 

(5^PA,5^P.)(fc•^) = iq^{z),kpx{z),r{z),kMz))- (249) 

The diagram 

TT*0 > T*Q (250) 



Tttq 
TQ 



TTQ 



is a vector fibration morphism. The vector fibration Tttq: TT*(5 — > TQ is dual to the vector fibration 
Ttq: TTQ TQ. The pairing 

(,)'^:TT*Q X TTQ^M (251) 

(Tttq, Ttq) 

is defined by 

{z,wr = ^{as),as))\s=o, (252) 

where C: R — > T*(3 and ^: M ^ TQ arc curves such that z = tC(0), w = t^(0) and ttqo( = tqo^. Such 
curves are provided by the coordinate constructions 

(9",Pa) o C = {q'^iz) + r{z)s,px{z)+px{z)s) (253) 

and 

(9«, q^)o^= {q'^iw) + q"'{w)s, q^{w) + q'^{w)s). (254) 
The coordinate expression of the pairing is 

{z,wr =p.{z)q'^{w) +Mz)rH- (255) 

A mapping 

V'Q:TT*g X T*Q^TT*g (256) 
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is defined by 



i^QiwJ) = w-tri{0) 



(257) 



with 77: R T*Q defined by r]{s) = rj*Q(ui) + sf. The coordinate expression of the mapping in terms 
of coordinates {q'^,Px,q^,'Pv) in TT*Q and coordinates {q'^, fx) in T*Q is given by 



{q'',p\,q^,p,^)o-tpQ = {q'',px,q'',Pu- f^)- 



(258) 



In the cotangent bundle T*TQ we use coordinates 



(g«,q\a^,6,):T*TQ 



J 4m 



(259) 



induced by coordinates {q'^,q ) in TQ. The LiouviUe form is the 1-form 



i^TQ = a^dq'^ + b^dq'^. 



(260) 



The 2-form 



t^TQ = da„ A dg** + d6« A d.j'' 



is the symplectic form on T*TQ. A vector fibration isomorphism 



(261) 



TT*g 

Tttq 
TQ 



= TQ 



(262) 



is defined as dual to the vector fibration isomorphism 



Kq 

TTQ ^ TTQ 



TQ: 



'TTQ 

= TQ 



(263) 



in the sense that 



{aQ{z),w) = {z,KQ{w)y 



(264) 



for z G TT*Q and w G TTQ such that T7rQ(2;) = tjq{w). We have the coordinate characterization 



{q'^A^,a^,b^) oaQ = {q'^,q^,p^,p„) 



(265) 



of the mapping aq. 
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For a vector z = tt^'^'^^x(0: 0) G TTT*Q represented by a mapping x- T*Q we have 

{aQl^TQ,z) = (?9tq,Tq;q(z)) 

= (^T*TQ(TaQ(2:)),T7rTQ(TaQ(z))) 
= (aQ(T-TT*Q(^))'TT7rQ(z)) 
= (^TT*QW''«QC'"T7rQ(z)))^ 

= (rT-r*Q(tt("'i)x(0,0)),«Q(TT7rQ(tt(°'i)x(0,0))))- 
= (t(0'i)x(0,0),«Q(tt("'i)(7rQox)(0,0)))~ 
= (tx(0,-)(0),tt(i'0)(7rQox)(0,0))- 

= ^(x(0,t),t(i-0)(^QOx)(0,i))|,.o 



dr 

_ _d 

~ dt' 

_ d 

~ di' 

d_ 
dt' 

(dT^?Q,tt(O'l)x(0,0)) 

{dT^Q,z). (266) 



-(x(0,i),t(7rQox)(-,i)(0))|t=o 
(^T*Q(tx(-, i)(0)), T7rQ(tx(-, i)(0)))|t=o 



= ;^T(^Q,t(^'°)x(0,t))|t=o 



We have used the formula (232) and relations 

"^T* TQ ° = aQO r-^T* q (267) 

and 

Ttttq o Tqq = TTttq (268) 

derived from 

TTTQ o ag = Tttq. (269) 

We have shown that 

a*Q^TQ = dT^Q- (270) 

It follows that the 2-form driOQ is non degenerate and that the mapping agiTT^Q T*TQ is a 
symplectomorphism from (TT*(5,dTWQ) to {T*TQ,cotq) since 

druiQ = drddq = ddr'&g = dag^jQ = agdt^jQ = Q:qWtq- (271) 
These results are confirmed by the coordinate calculations 

a*Q^TQ = P^dg" + p^dr = dr^g (272) 

and 

agUTQ = dp„ A dcf + dp„ A dg" = dr^g. (273) 
In the cotangent bundle T*T*Q we use coordinates 

(g«,PA,n^,t;^):T*T*Q^K4'" (274) 

25 



induced by coordinates {q'^,p\) in T*Q. We have the Liouville form 



and the symplectic 2-form 



0Jj*Q = duK, A dg" + dv'^ A dp^,. 



on T*T*Q. We have aheady introduced the mapping 



characterized by the equahty 



(/3(T*Q,a,Q)W'^) = (WQ,UAt;) 

for vectors u G TT*Q and v G TT*Q such that Tj*q(u) = Tj*q(u). The diagram 



^(T*Q,a.g) ^^^^^ 



= T*Q 



TT*Q 

is a vector fibration isomorphism. For each z G TT*T*Q we have 

= (^T*T*Q(T/3(T*Q,c.Q)(^))'T7rT*Q(T/3(-r*Q_^^-,(z))) 

= (/3(T*Q,^Q)('rTT*QW)'TrT*Q(z)) 

= ('^Q.Ttt*q(^) Ttt-*q(2)) 

The formula (218) and relations 

Tt*T*Q ° "''/?(T*Q,t.Q) = /5(T*Q,a;g) ° ^TT*Q 

and 



derived from 



were used. We have shown that 



(275) 
(276) 

(277) 
(278) 

(279) 



(280) 

(281) 

(282) 
(283) 
(284) 



It follows that the mapping /3(j*q TT Q T T Q is a symplectomorphism from (TT Q, dpWQ) 
to (T*T*<5, Wy*g) since 

drWQ = diTWQ = d/3jV*Q_^^^ Vq = /^(T*Q,c.o)dVQ = /5(T*Q,a,g)'^T*Q- (^^S) 

Coordinate calculations 
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(286) 



and 

= dp^ A dg« + dp^ A dr = druJQ. (287) 

confirm these results. 

7. Dynamics of mechanical systems. 

Let Q he the configuration manifold of a mechanical system. The cotangent bundle T*Q is the phase 
space of the system. Elements of the phase space are momenta. The commutative diagram 




(288) 



Q 

contains the geometric structures used to formulate the dynamics of the system. The dynamics is a 
differential equation D c TT*Q. A solution tt: 7 ^ T*Q of this equation is a phase space trajectory of 
the system. External forces have to be included in a complete description of dynamics. The dynamics 
of the system with external forces is the differential equation 

Df =tljQ\D) CTT*Q X T*Q. (289) 

A solution is a curve (tt, (p): I T*Q x J*Q. The values of this curve represent the momenta and 

external forces. The differential equation is of first order for the momentum component tt and of 
zero order for the force component tp. This treatment of external forces is suitable for non relativistic 
systems. Relativistic systems described by homogeneous Lagrangians require a modification of the 
concept of force. We will deal with dynamics without external forces. 

Trajectories of the system in the configuration manifold Q are solutions of the second order Euler- 
Lagrange equation 

E = TVq(P£>). (290) 

We have recognized the presence of a canonical symplectic structure in TT*Q with the symplectic 
form dxoJQ ■ In most cases of interest in relativistic physics the dynamics is a Lagrangian submanifold of 
{TT*Q, dr^g). Morphisms uq and /3(t*q are canonical symplectomorphisms from (TT*Q, drwg) 
to {T*TQ,ujtq) and to (T*T*(5, w-|-*q). These symplectomorphisms with cotangent bundles create 
the possibility of generating the dynamics from (generalized) Lagrangians associated with TQ or (gen- 
eralized) Hamiltonians associated with T*Q. 

We will present a number of examples of mechanical systems in Lagrangian and Hamiltonian for- 
mulations. We will perform the Legendre transformations and test the integrability criteria for these 
systems. 

8. Lagrangian systems. 
Let 

L:TQ^R (291) 

be the Lagrangian of a mechanical system with configuration space Q. The Lagrangian may be defined 
on all of TQ or on an open subset. The image N = im(dL) of the mapping 

dL: TQ T*TQ (292) 
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is a Lagrangian submanifold of {T*TQ,ujtq) and the set D = aQ^{N) c TT*Q is a Lagrangian 
submanifold of {TT*Q,dTU)Q). Coordinates {q'^,q^,an,bi,) of N satisfy equations 



b^ = di,L{q\q^) 



(293) 



and coordinates {q'^,Px,q^,'Pv) of elements of D satisfy equations 



(294) 



p^ = a^L(g^g^) 
derived from the variational principle 

pjq'' +pjr = SL{q-,q^) = d^L{q-,q^)Sq>' + d(,L{q'^,q^)5r. (295) 
Substituting the equalities 



and 



dp^ = d^d^Ldq" + dcd^Ldr (296) 

dp^ = d^df^Ldq" + di,df,Ldr (297) 

in 

drcoQ = dpi, A dq" + dp,, A dq" (298) 

we obtain the equality 

dTLJQ\D = d^d/^Ldq" A dq^ + di,d^Ldq" A dqi^ + d^df.Ldq" A d^** + dodf.Ldq" A dg" = 0. (299) 

This equality together with dim(_D) = 2m = 1/2 dirn(TT*Q) c;onfirrns that _D is a Lagrangian subman- 
ifold of {TT*Q, dxcOQ). The set is a Lagrangian submanifold since it is the image of the differential 
of a function and D is a Lagrangian submanifold since it is obtained from N by applying the symplec- 
tomorphism ag. We have confirmed this fact by direct calculation. 

The set D C TT*Q is a differential equation. A solution is a curve 7:M — » T*Q such that vectors 
tangent to 7 are in D. Equations 

7. = a,L(7^7^) 
are differential equations for the coordinate expression 

(7",7a) = (9",Pa)o7 (301) 

of a curve 7 derived from equations (294). Dots indicate derivatives. The differential equation D 
represents dynamics in the sense that solution curves are phase space trajectories of the mechanical 
system. We say that Z) is a Lagrangian system since it was obtained from a Lagrangian function (291). 
Equations (294) and equations (300) will be called the Lagrange equations. 
The second order equations 

p^ = d^L{q-,q^) 
p^ = df,L{q-A^) 

P^ = a.9^L(g^ q^yt + d„d^L{q-, q^W 

p„ = d,df,L{q-, q^W + dodf,L{q\q^W 
represent the prolongation PD of the Lagrange equations. The equations 

a,5^L(g^ + d,df,L{q\ q^W - d^L{q-, q^) = (303) 
are the Euler-Lagrange equation T^7rQ(P£>) in coordinate form. 
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(302) 



Example 4. Let Q be a manifold of dimension 3 with coordinates (g*) = {q^,q^,q'^)- There is a 
Riemannian metric tensor gij, a function (p, and a 1-form A = Ajdg' on Q. The function 

L{q\ q') = -^9ij€Q' -e^ + eAiH (304) 
is the Lagrangian of a particle of mass m and charge e in an electric field 

E = Eidq' = -dip = -dupAq^ (305) 

and a magnetic field (induction) 

B = iSydg* A dqi =dA = dA, A dq' = d^Ajdq' A dq^ = ^{diAj - djAi)dq' A dqK (306) 
The dynamics D of the particle is described by equations 

Pi = ^digjuifct - ediip + ediAjq' ^^^^^ 

Pi = mgijq^ + eA^ 

Gauge independent covariant second order Euler-Lagrangc differential equations 

mgijiij^ + V^q'^q^) = eEi + Bijq^ (308) 
are easily derived. The symbol T^^ is the Christoff'el symbol 

^ki = Id'' {dkgii + dighi - digki ) • (309) 

Solution curves of (308) are motions in the configuration space Q. These equations provide a partial 

description of dynamics. The complete description of dynamics is obtained by complementing these 
equations with the gauge dependent velocity-momentum relation 

Pi = mgijq^ + eAi. (310) 

A gauge transformation 

Ai^Ai + diip (311) 

will not modify equations (308) but will change the velocity-momentum relation. ▲ 

Example 5. A gauge independent formulation of dynamics of a charged particle is obtained by 
extending the configuration space Q to a manifold Q of four dimensions with coordinates {q,q^). A 
gauge transformation is a coordinate transformation (g, q^) ^ {q-\-'4^{q^), 3*)- The dynamics is derived 
from the gauge independent Lagrangian 

L{q, q\ q,q') = -^gaq'q' - et^ + eAiq' + eq. (312) 

The equations 

p = 

Pi = ^digjkq'q^ - edup + ediAjq> 

^ (313) 
p = e 

Pi = mgijcf + eAi 

provide a description of dynamics in terms of coordinates {q,q^,p,pj,q,q'^ ,p,pi) in TT*Q. These 
equations are gauge independent and can be given an explicitly covariant and gauge independent form 

p = 

mgij {q' + Tliq'^q') = eEi + Bijq^ ^^^^^ 
p = e 
Pi - eAi = mgijq>. 

The gauge invariant quantity [pi — eAi) is the momentum of the particle. ▲ 
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Example 6. Let Q be the space-time of general relativity with coordinates (g*^) = g^, g^, g^). 
The gravitational field is represented by a Minkowski metric g^x and the electromagnetic field is a 
2-form 

F = ^F^xdq'' A dg^ = -dA = -dAx A dg^ = -d^A^dq'^ A dg^ = ^{dxA, - d^Ax)dq'' A dg^ (315) 
derived from a potential 

A = A^dq"-. (316) 
The dynamics of a relativistic particle of mass m and charge e is derived from the Lagrangian 



L(g^ g^) = m^/g^xq'^q^ + eA^q'^ (317) 

defined for time-like vectors - vectors satisfying QKxq'^q^ > 0. Dynamics D C TT*(5 is described by 
the Lagrange equations 

= md^gxA;^^ + eS^^Ag^ = mgxJ^^AAr + ed^Axq^ 

m\\ m (318) 

Pk = mg^x— + eA^ 

\m 



with ||g]| = ^/OkxA^P^- Note that these equations are rcparamctrization invariant: if 7:R ^ T*Q is a 
solution and u: R ^ ffi. is a diffeomorphism with positive derivative, then 7 o a is a solution. One can 
say that solutions are one- dimensional oriented but not parametrized submanifolds of T*Q. 
The Euler-Lagrange equations 

^5«A (9^ + ^^.rr) + eF^^e = jr^g.xq^g^. + TpJ^^) (319) 

are reparametrization invariant and gauge invariant. If proper time is chosen as the parameter, then 
||g|| = 1 and the world line of the particle in space- time is a solution of simplified gauge invariant 
second order differential equations 

m{r + r^^q'^r) = -eg^^Fx^r- (320) 

The complete dynamics is gauge dependent. A 

Example 7. Let Q be the space-time of general relativity with coordinates (g'' ), a Minkowski metric 
g^x and an electromagnetic potential A = A^dq'^ . Let Q be a manifold of dimension 5 with coordinates 
(g, g""). A gauge transformation is a coordinate transformation {q,q'^) 1-^ (g + i^{q^),q'^)- We use 
coordinates {q,q'^,p,Px) in T*(5. Two gauge invariant quantities arc derived from the ^-momentum 
{p,px)- These are the charge p and the energy momentum (p\ — pA\). There are two equivalent 
interpretations of the manifold Q [6]. This manifold is interpreted as a pseudoriemannian manifold 
(Kaluza [5]) with a metric tensor 

■ 1 Ax 

A^ g«A + A^Ax 



(321) 



or as the total space of a principal fibration (Utiyama [14]) 

X:Q^g (322) 

characterized by 

{q'')ox={q^) (323) 
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The electromagnetic potential is used to introduce the connection form 

a = dq + A^dq" (324) 
in the principal bundle Q. The curvature form 

f3 = da = -^F^xdq'^dq^ = dAx A dq^ = d^Axdq'' A dg^ = -^{dxA, - d^Ax)dq'' A dg^ (325) 

represents the electromagnetic field. 

The Lagrangian of a particle with mass m and charge e is the gauge invariant function 



(327) 



L{q, q'^, q, q^) = my/g^xq^q^ + eA^q" + eq. (326) 
Coordinates {q,q'^,p,Px,q,q^,P,'Pv) are used in TT*Q. The Lagrange equations 

p = 

p„ = "^^«fi'AA.||||||- + ed^Axq'' = mgxuTl^i^^^ + ed^Axq^ 
p = e 
= mg^x-r-TTT + eA^ 

\m 

are equivalent to the explicitly covariant and gauge independent second order equations 

p = 

m{r + r^^q^qn = eg^^Fx^q^ 

p = e 
Pk - eA^ = mg^xq^- 



(328) 



These equations arc obtained by adopting the simplifying condition ||g|| = 1. Trajectories in Q satisfy 
the second order equations 

m(r + Vx^q^r) = eg'^^Fx^r (329) 

with no conditions on q. Compatibility with field equations requires that trajectories in Q be two 
dimensional. The dynamics of a charge particle has to be suitably modified for correct description of 
interaction with the electromagnetic field. ▲ 

Not all mechanical systems are Lagrangian systems derived from a Lagrangian defined on the tangent 
bundle TQ. The dynamics could be generated by a Morse family of functions defined on fibres of a 
fibration 

r]:Y ^ TQ. (330) 
We have coordinates {q'^,q^) in TQ. In the space Y we use adapted coordinates 

(g'^,g^,y^):F^M2m+fc (331) 

such that 

(9^9^)or/ = (g^g^). (332) 
Let L: F ^ R be a Morse family of functions defined on fibres of r]. The k x (2m + k) matrix 



/ d^L d^L d'^L 

\ dy^dy^ dy^dq'^ dy^dq^ 



(333) 
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is of maximal rank. The Lagrangian submanifold N c T*TQ generated by the family is the set of 
elements of T*T(5 with coordinates [q"' ,q^,an,hv) satisfying equations 

a^ = d^L{q'^,q\y^) 

K = di,L{q\q^,y^) (334) 
Q = dAL{q\q^,y% 

for some values of the variables y^. The set D = a^{N) C TT*Q is a Lagrangian submanifold of 
(TT*(5, dr^g). It is the set of elements of TT*(5 with coordinates {q'^,px, q^,Pu) satisfying equations 

= 9,L(g^g^^/^) (335) 
= dAL{q^,q\y^) 

for some values of the variables y^. These equations can be derived from a variational principle 

P.Sq^+P.Sq- = SL{q-, q\ y^) = d^L{q-, q\ y^)dq^+df,L{q-, q\ y^)Sr+dBL{q-, q\y^)Sy''. (336) 

Equations (335) present the set D parametrized by variables {q'^, q^, y^) subject to dAL{q'^, ,y'^) = 0. 
From 

dp^ = d^d^LAq'' + d^d^LAq'' + dAd^LAy^, (337) 

dp^ = d,df,L&q'' + dc.d(,Ldr + dAd(,Ldy^, (338) 

OaO^L = 0, (339) 

dAdf^L = (340) 

we obtain the equality 

dTiOQ\D = 0. (341) 

Equations dAL{q'^ ,q^ ,y^) = leave only 2m out of the 2m + k variables {q'^,q^,y"^) independent. 
This is a consequence of maximality of the rank of the matrix (333). It follows that dim(D) = 2m. We 
have thus confirmed that £> is a Lagrangian submanifold of (TT*Q, drwg). The set £> is a differential 
equation and may represent the dynamics of a mechanical system. 

Example 8. Let Q be the space-time of general relativity with coordinates {q'^) = (g°, q^,q^,q^) and 
a Minkowski metric g^x- Let 

L{q\q\y) = ^9.xq''q^ (342) 

o o 

be a function on TQ x K_|_ , where TQ is the tangent bundle with the zero vectors removed. This function 
is a Morse family of functions of the variable y with the coordinates [q"" , q^) treated as parameters. It 
represents the Lagrangian of a particle of mass zero. The dynamics of the particle is governed by the 
equations 

Pk = TT^kgxnq^q^ 
2y 

Pk = '-9Kxe ^'''^ 

y 

= gKXq'^q^ 

satisfied for some value of the variable y. The variable y can be eliminated from the equation for p^. 
It follows from the resulting equation 

Pk - T^^pxr = (344) 
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that the covector p„ is covariant constant along the world line. If an afRne parameter is chosen then 
y is constant and the dynamics satisfies equations 

(345) 



1 X 

Pk = -g^xq 

y 

for some constant y > 0- A 

9. Hamiltonian systems. 

Let 

H:T*Q^R (346) 
be the Hamiltonian of a mechanical system with configuration space Q. The mapping 

-dH:J*Q ^J*J*Q (347) 
is a section of the fibration tTj^q. Consequently 

^ = ^(T*Q..p) ° (-d^): T*Q ^ TT*Q (348) 

is a vector field. The image M = im(— diJ) is a Lagrangian submanifold of {T*T*Q,u!j^q) and 
D = im{X) is a Lagrangian submanifold of (TT* Q , (Ituiq) . The equations describing D are the 
Hamilton equations 

q'^ = d'^H , ^ 

derived from the variational principle 

p^Sq^^ - rSp^ = SH{q-,px) ^ -d^H(q-,px)Sqf^ - d^H(q\px)Sp^. (350) 
The symbol d'^ denotes the partial derivative . The dimension of D is 2m and by using equalities 

dq" = d^d^Hdq" + d^d^Hdp^ (351) 

and 

dp^ = -d^d^Hdq" - d'df.Hdp, (352) 

in 

druJQ = dp^ A dq^ + dp^ A dq^" (353) 

we obtain the equality 

dTUjQ\D = -d^d^.Hdq" A dq" - d'^d^.Hdp^ A dq" + d^d^Hdp^, A dg^ + d^'d^Hdp^, A dp^ = 0. (354) 

It follows that D is a Lagrangian submanifold of (TT*(3, drwg). The set is a differential equation 
and may represent the dynamics of a mechanical system. 

Example 9. Equations (307) of Example 4 can be rewritten in the form 

q' = ^g'HPj - eAj) 

I g (355) 

Pi = ^digjug^^g'^'^iPi - eAi){pm - eA^) - ediip + —diAjg^\pi - eAi). 

These equations describe a Hamiltonian vector field. They are the Hamilton equations for the Hamil- 
tonian 

H{q\pj) = ^g''{Pi - eAi){pj - eAj) + e^. (356) 
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Gauge independent dynamics of charged particles and the dynamics of relativistic particles are not 
images of Hamiltonian vector fields. Dirac [1] introduced generalized Hamiltonian system,s in order 
to be able to deal with similar cases. In the original construction of Dirac a generalized Hamiltonian 
system is a family of Hamiltonian vector fields on the phase space T^Q restricted to a constraint set 
C C T*Q. We have translated this construction in an equivalent construction of a differential equation 
D c TT*Q. 

Let C C T*Q be a submanifold and let 

H-.C^m. (357) 

be a differentiable function. The set 

M={bG T*T*Q; a = n^^Q{b) e C, V„^-r<.ccT„T*Q<^' = -i^H,u)} (358) 
is a Lagrangian submanifold of (T*T*(5, Wj*g) and 
D = 0-^\^^^JM) = {we TT*Q; a = Tj*q{w) e C,\/^^^^^^^^^^^{ujq,uAw) = {dH,u)} (359) 

is a Lagrangian submanifold of (TT*(5, d^WQ). If is a set of k independent functions on T*Q 

such that 

C = {ae T*Q; $^(0) = for ^ = 1, . . . , fc} (360) 

and i? is a fimction on T*Q such that H\C = H, then coordinates {q'^,P\,q^,Pi^) of elements of D 
satisfy the equations 

(f = d^H + v^d^^A (361) 
p« = -dji + v^d^^A 

derived from the variational principle 

p,S<f-rSPf = SH{<f,P>.} = -%H{g',px)Sqi'-a'H(q',p^)Sp^ 
with variations {6q'^,5px) satisfying 

5«$^(5,7« + d''<i>ASp^ = 0. (363) 

Lagrange multipliers {v^ G M.'^) appear in these equations. The same equations are derived from the 
variational principle 

P^.Sq'' - rSpf. = SH(q\px,v'') = -df,H{q'',px)6q'' - d^H{q\px)6p^ + ^ Aiq"" ,Px)5v^ (364) 
corresponding to the Morse family 

H{q'^,px,v^)=H{q^,px)-v^^A{q^,Px) (365) 
of functions of the variables [v^ e K*^). The set D is & Lagrangian submanifold since it is generated 



by a Morse family. At each point a G C the set Da = D (1 TaT*Q is an afRne subspace of the vector 
space TaT*Q. 

In Dirac's construction the dynamics is described by vector fields 

X = o {v^d<^A -dH):T*Q^ TT*Q (366) 

with arbitrary functions v^{q'^,px)- These fields are restricted to the constraint set C. 
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Example 10. The dynamics of the charged particle of Example 5 is a Dirac system. The Hamiltonian 
is the function 

H{q, q\p,Pj) = ^9'^(Pi - eAi){pj - eAj) + etp (367) 

and the constraint is the set characterized by 

^q,q\p,Pj)=p-e = 0. (368) 

The function ^ 

H{q, q\p,Pj, v) = —g'^ {pi - eAi){pj - eAj) +eip- v{p - e) (369) 

is a Morse family of functions of e M. Equations 

p = e 
q = V 

= - eAj) (370) 

p = 

1 e 
Pi = T^digjkg^^g''"'(pi - eAi){pm - eA„) - edup + —diAjg^\pi - eAi). 

2 777, m 

obtained with this Morso family arc equivalent to equations (313). A 

The dynamics of a non relativistic charged particle in the above example is the only Dirac system 
known to us. Hamiltonian formulations of relativistic dynamics require a higher level of complexity. 
Differential equations generated by Lagrangians and by Lagrangian Morse families are Lagrangian 
submanifolds of (TT*(5, dpWQ). The same is true of Dirac systems. We define a generalized Dirac 
system as a differential equation D c TT*Q, which is a Lagrangian submanifold of {T'V*Q,At<jJq)- 
Existence of Morse families for open subsets of Lagrangian submanifolds is guaranteed by Hormander's 
theorem. Known generalized Dirac systems are globally generated by Hamiltonian Morse families. 

Example 11. Lagrange equations (318) of Example 6 have an equivalent form 
g'^^iPK - eA^){px - eAx) = rr? 

= - (371) 

V V€ 

Pk = -7^d^g^^{p^ - eA^){p, - eA,) + —g'"'d^A^{p, - eA,) 
2m ft- m 

with arbitrary v > Q. These equations are obtained from the Morse family 



H{q^,px,v) = V \y^Jg''^{pn - eA^)(px - eAx) - mj (372) 

of functions of the variable v > 0. A 

Example 12. The gauge independent dynamics in Example 7 is a generalized Dirac system. Equa- 
tions 

g'^^ip^ - eA^){px -eAx) = vr? 

p = e 
q = v^ 

r = -9^\PX - eAx) (3^^) 
m 

p = 

1 1 

Pk = -^d^g^'-'ip^ - eA^){p^ - eA^) + —g^'''^^A^,{Jp^ - eA^) 
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with > and arbitrary v"^ are equivalent to equations (328). These equations are generated by the 
Morse family 



H{q'',px,v\v^)=v^ J g'^^ip^ - eA^){px - eAx) - m] +v^p-e) (374) 



of functions of the variables > and w^. ▲ 
Example 13. The generalized Dirac system of Example 8 is described by equations 

5" Vpa = 

= vg^^px (375) 

V 

derived from the Morse family 

H{q'',px,v) = ^g'^^p,px (376) 
with V > 0. A 

10. The Legendre treinsformation. 

A Lagrangian L{q'^, q^) is said to be hyperregular if the Legendre mapping 

A:TQ^T*Q (377) 

defined by 

{q-,px)oX = {q-,d^Liq'',r)) (378) 
is a diffeomorphism. Let the mapping 

X:T*g^Tg (379) 

represented by 

(<^^g^)ox = (g^x^9^Pa)) (380) 

be the inverse diffeomorphism. Relations 

df,L{q^,X-{q',Pa))=P. (381) 

x'^(g^9^'L(g^«'^))=r (382) 
hold. Using the diffeomorphism x to eliminate the velocities {q'^) from the energy function 

E{q^,Px,r)=Pxq^-Liq'',r) (383) 

we obtain the Hamiltonian 

H{q^,p.)=Pxx\q^,Pu) - Liq^,X-'iq',Pa)). (384) 

The energy function is defined on T*Q x TQ. The passage from the Lagrangian to the above 

Hamiltonian is the Legendre transformation for a hyperregular Lagrangian. 

Let H be the Hamiltonian (384) obtained from a hyperregular Lagrangian L. For the mappings 

d7f:T*Q^T*T*Q, (385) 
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and 

V*Q,.,)°(-d^f)^T*Q^TT*Q (387) 

we have 

{q\px,u^,vnodH =iq-,px,d^HiqP,p„),d-'HiqP,p,)), (388) 
(<?^PA,g^P.) o/3-^i^^_^^j = {q\px,-v'',u.), (389) 

and 

{q^Px,q^,Pu)of]-^\^^^^^ o i-dH) = {q-,px,d^'H{q'>,p,),-d,H{qP,p„)). (390) 
On the other hand, we have 

(g«, q\ a^, K) o dL = (g«, g\ a^i(g^ g^), di,L{qP, ?")), (391) 
((Z",PA, 9",^.) o ag' = (?", r, a.), (392) 
(9",PA, ?^P.) o o di = (g«, a^L(9'', g^), g^ a,L(g^ ^ )), (393) 

and 

(g^PA,4^P.)oaQlodLox = (g^5J,i(g^x"(g'^,P.)),x''(a",Px),5.L(g^x"(a",Pr))) (394) 
for the mappings 

dL:TQ^T*TQ, (395) 
ag^T*TQ^TT*g, (396) 
ag^ o dL: TQ ^ TT*Q, (397) 

and 

ag^ odLox:T*g ^ TT*Q. (398) 

From 

a«^^(g^p,) = p^a«x''(9^P.) - 5«i(g^x''(g^p.)) - api(g^x"(9^pr))a.x''(g^p.) 

= -a«i(g^x"(9^Pa)) (399) 

and 

d^H{q>^,p,) = x''{q',P.)+P^.d^X^{qP,p,) - dpL{q^ ,x%q^ ,Pr))d^X%q'' ,Pa) 

= x"(a^Pa) (400) 

it follows that 

V°diox = /3-^.g_^^)0(-diJ) (401) 
We see that the Hamiltonian and the Lagrangian generate the same dynamics 

D = MP^j\q^^^^ o i-dH)) = im(agi o dL). (402) 
Example 14. The Lagrangian 

77% 

L{q\ q') = -^Oijq'q' - + eA^q' (403) 
of Example 4 is hyperregular. The Legendre mapping and its inverse are the mappings 

{q\Pj) o A = {q\mgjkq'' + eAj) (404) 

and 

{q\ q')ox= {q\ ^g'' {Pj - eAj)j . (405) 



From the energy function 



E{q\pj, q'') = pd - jQijq'q^ + ey. - eAiq' (406) 



we derive the Hamiltonian 

H{q\pj) = ^g'^iPi - eAi)ipj - eAj) + ecp. (407) 



▲ 
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If a Lagrangian is not hyperregular, then the Legendre mapping is not invertible. It may happen 
that the image of the Legendre mapping A is a submanifold C C T^Q. Let G (lT*Q x TQ be 

the graph of A (intersected with T*(5 x TQ c T*(5 x TQ). The energy function (383) restricted 
to the graph G does not depend on q'^ since 

dkE =p^- df,L. (408) 

If fibres of the projection prf.T^Q x TQ ^ T*(3 are connected, then a function 

H{q-,px) = E{q\p^,r) (409) 

can be dcfin(;cl on C by substituting in E any values of the velocities (g^) such that — d^L = 0. 
This is the construction of the constrained Hamiltonian used by Dirac. This is also the construction 
of the generalized Legendre transformation introduced by Cendra, Holm, Hoyle and Marsden. The 
Dirac system generated by the constrained Hamiltonian is sometimes equal to the Lagrangian system 
generated by the original Lagrangian. 

Example 15. Applying the Cendra-Holm-Hoyle-Marsden generalized Legendre transformation to 
the Lagrangian 

Tfi 

q\ q, ct') = jgzjg'g^ - e<^ + eAiq' + eq. (410) 
of Example 5 we obtain the constraint C described by 

p = e (411) 

and the Hamiltonian 

H{q, q\Pj) = ^g'HPi - eAi){pj - eAj) + e^p (412) 

defined on the constraint C with coordinates {q, q'',pj). This Hamiltonian generates the equations (370) 
of Example 10 equivalent to the Lagrange equations (313) of Example 5. The Cendra-Holm-Hoyle- 
Marsden version of the Legendre transformation gives correct results for this example of a mechanical 
system. ▲ 

The mechanical system in the above example is the only mechanical system known to us for which 
this version of the Legendre transformation functions correctly. 

Example 16. The Lagrangian 



Liq", q^) = mVg^oA^ + eA^q" (413) 
of Example 6 is singular. The image of the Legendre mapping 

{q'',Px)o\={q'",mgx^,-^^+eAx) (414) 

is the constraint set C described by 

g''^{p^-eA^){px-eAx)=m'. (415) 

The energy function 

^(9%PA, q") = Pxq^ - m^/g^x^q^ - eA^^ (416) 
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vanishes on the graph G of the Legendre mapping and the Dirac Hamiltonian is zero. Differential 
equations derived from this constrained Hamiltonian are the equations 



9'^^{Pk - eA^){px - eAx) = rr? 

r = ^5"^(PA - eAA) (417) 



m' 



Pk = -T^d^g^'^ipi, - eA^){p, - eA,) + —g^^'^d^A^{p, - eA,) 

with arbitrary values of the Lagrange multiplier v. The Dirac system represented by these eqiiations 
is the union Ufg UD- of three sets corresponding to > 0, = 0, and v < respectively. The set 
Z)+ is the generalized Dirac system D of Example 11 equivalent to the Lagrangian system of Example 
6. The set Do described by 

g^'^iPK - eA^){px - eAx) = vr? 

r = (418) 

must be excluded since the velocities qn evaluated on vectors tangent to world lines are never zero. 
The set Z)_ is a generalized Dirac system obtained from the Lagrangian 

L.^q^'A^) = -m^g^xq^q^ + eA^q^. (419) 
Setting 11 = 1 in equations (417) we obtain equations 

Pfi = mg^xq^ + eA^ 

(420) 

m(r + r^^q'^r) = -eg^^Fx^r 

correctly describing the dynamics of charged particles. Solutions of these equations are curves using 
proper time as the parameter. With y = —1 equations (417) result in the equation 

p^ = -mg^xq^ + eA^ (421) 

and the second order system 

m(r + T^^q^q^) = eg^^Fx^q^. (422) 

Solutions of the second order equations are world lines of particles with mass m and charge — e or 
particles with mass — m and charge e. The equation (421) suggests that we are dealing with particles 
with negative mass — m. The principle that world lines of particles with positive energy should be 
oriented towards the future and that world lines of particles with negative energy (antipartieles) should 
be oriented towards the past (Stueckelberg [11], Feynman [3]) is violated. We conclude that the Cendra- 
Holm-Hoyle-Marsden version of the Legendre transformation is too fast to provide correct results for 
this important example of a mechanical system. ▲ 

The Dirac system generated by a constrained Hamiltonian 

H:C ^M. (423) 

is characterized by the variational relation 

pjq'' - q^5p^ = -d.kHSq'' - d^H6p^ (424) 

on C. A constrained Hamiltonian derived from a singular Lagrangian can be considered a function on 
the graph G of the Legendre mapping defined by 

H{q'',px) = E{q\px, q^ = px^" - L{q\ r) (425) 
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The variational relation 



pM"" - q^Sp^ = -d,kE6q^ - d'^ESp^ - d^Edq'' (426) 

on G is equivalent to the relation (424) for the Hamiltonian (425). The variations {6q'^,6px,Sq^) are 
components of a vector tangent to G. Hence 



If Lagrange equations 

are satisfied, then 
and 



Sp. = d^dk6q>' + df,dk6r. (427) 

p^ = d^L{q-,q^) 



X, (428) 



p, = 5^L(g^g^) 

pjq'' - rSp. = d^LSq'' - ^Sp. (429) 



-d,kE6q^ - d^ESp^ - dkESr = d^LSq'^ - ^Sp^ - (p« - dkL)Sr = d^L6q^ - q^Sp^- (430) 

This seems to imply that the Dirac system generated by the Hamiltonian (425) is equivalent to the 
Lagrangian system (428). This conclusion is not correct. It is true that the variational relation (424) 
follows from the; Lagrange equations. The converse is not true since the same constrained Hamiltonian 
may be in the relation (425) with different energy functions constructed from different Lagrangians. 
The Lagrangians L+ = L and L_ of Example i6 generate different Lagrangian systems and £)_ 
but lead to the same constrained Hamiltonian. 

We observe that if E(q'^,p\,q'^) is the energy function associated with a Lagrangian L(q'^,q'^), then 
E[q'^ ,p\,v^) is a Morse family of the variables (u^). The rank of the matrix 



d'^E d'^E d^E \ f ^d^L -d'^L 



dv'^dv^ dv'^dqt^ dv'^dpi, J \ dv'^dv^ dv'^dqi^ 



(431) 



is maximal. No requirements are imposed on the Lagrangian. The generalized Dirac system generated 
by the Morse family E{q'^,p\,v^) is obtained from the variational relation 

pjq^ - r^PK = -d,kE6q'' - d^E5p^ + df^ESv^ (432) 

with arbitrary variations {6q'^ , Spx, Sv^) . With 

E{q-,px,v^^)=Pxv^-L{q-,v^) (433) 

the relation takes the form 

pjq- - rSp. = d,L{q\ v^^)Sq'' - v'^dp, - {p, - 5^L(g^ v^^))6v\ (434) 

Equations 

p^ = d^L(cf,v^) 

p, = 5^L(g^^;^) (435) 
q^ = v'^ 

obtained from this relation are equivalent to the Lagrange equations. 

The passage from the Lagrangian L{q'^,q^) to the global Hamiltonian Morse family E{q'^ ,px,v^) is 
the slow and careful Legendre transformation required to provide a correct Hamiltonian formulation 
of any Lagrangian system. The Morse family generating a Lagrangian submanifold is not unique. 
Modifications resulting in a reduction of the number of variables are usually possible. 
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Example 17. The Hamiltonian Morse family 

E{q,q\p,Pj,v,v^) =pv+piv'' - ^gijV^v^ + eip - eAiv'' - ev (436) 
for the charged particle of Example 5 reduces to the Hamiltonian Morse family 

H{q, q\p,pj,v) = -^9"'{Pt - eAi){Pj - eA,) + eip ~ v{p - e) (437) 
of Example 10. The reduction is obtained by using in (436) the equality 

v"" = -g'^^Pj - eAj) (438) 
m 

obtained from 

dE 

— =Pi- mgijV^ - eAi = 0. (439) 

▲ 

Example 18. 

E{q'^,px, v^) = p^v" - m^/g^xV^v^^ - eA^v'^ (440) 
for the charged particle of Example 6 reduces to the Hamiltonian Morse family 

H{q'',px,v) = V (^Y^5«^(p« - eA^){px - eAx) - (441) 

of a single variable v > introduced in Example 11. Variables {q'^,px,v^) are coordinates in the 
manifold T*Q x JQ. Only the open subset described by the inequality g^xv^'v^ > is considered. 

Variables {q'^,Px,v) are coordinates in T*Q x M+. The function E{q'^,px,v^) is a Morse family of 
functions on fibres of the fibration 

C:T*Q X Tg^T*Q (442) 

characterized by 

{q^,Px)oC = {q\px). (443) 
This fibration can be interpreted as the fibration 

C':T*Q X Tg^T*gxM+ (444) 

characterized by 

{q\px,v) o C' - {q\px, ^yg^uv^'v-) (445) 

followed by the projection 

prT*g:T*QxR+^T*Q. (446) 
Fibres of (' are hyperboloids g^vV^v'^ = . From 

bE{<t,px,v>') = -^5v'' = ^ (447) 

with variations ^v^ satisfying 

Ka^j.^v^v") = Ig^.v^bv'' = (448) 

we obtain two critical points 

zS'^PA-e^A) (449) 



^g^^^vvi^v^- 

in a fibre over the point with coordinates {q'^,Px, v). Evaluating the function E at these points results 
in two Morse families 



H±{q^,px,v) = ±v \^.^g^^{p^-eA^){px-eAx)Tmj (450) 

of functions of v. One of these is the reduced Morse family (441). The other generates an empty set 
since it has no critical points. The family (441) depends linearly on the variable v restricted to positive 
values. No further reduction is possible. ▲ 
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Example 19. The dynamics of the charged particle of Example 7 is generated by the Hamiltonian 
Morse family 

E{q,q'^,p,P\,v,v'^) = pv + Pf.v'^ - g^xV^v^ - eA^v"^ - ev (451) 
or by the Hamiltonian Morse family 



H{q'',px,v^,v'^)=v^ (jg'^^Pn - eA^){px - eAx) - +v^{p-e) 

of functions of two variables > Q and v'^ as in Example 12. 

The slow Legendre transformation can be extended to Lagrangian Morse families. If 



is a Morse family, then the A; x (A; + 2m) matrix 



(. 



\ dy^dy^ dy^dv^ dy^dp^ J 
is of maximal rank k. It follows that the function 

E{q'^,px, t;") = p.v'^ - L{q'^, v\y^) 
is a Morse family of functions of the variables {y^, v'^) since the (A: + to) x [k + 3to) matrix 



(452) 



(453) 



(454) 



(455) 



d'^E d'^E 



dy^dy^ dy^dv^ dy^dq^^ dy^dpi, 
d^E d^E d'^E d^E 



dy^dy^ dy^dv^ dy^dqf^ 

-d'^L -d'^L -d^L 



dv'^dy^ dv'^dv^ dv'^dqi^ dv'^dpi, ' ^ dv'^dy^ dv'^dv^ dv'^dqi^ 
is of rank k + m. 

Example 20. The Hamiltonian Morse family 

E{q'',Px, y, w^) = p^v" - ^g^xv^v^ 

generates the generalized Dirac system of Example 8. The equations 

di,E{q'',px,y,v^)=Pi, - ^9n\v^ 

permits the elimination of v'^ from E{q'^,px, y, v^). The result is the Hamilton Morse family 

H{q^,px,y) = I^'Vpa 

with y > 0. It is the Morse family of Example 13 with v replaced by y. 

42 



(456) 



(457) 



(458) 



(459) 



11. Integrability. 

In Section 4 we have established intcgrabiUty criteria for a class of differential equations which can be 
specified in terms of vector fields restricted to a submanifold. The dynamics of relativistic mechanical 
systems presented in our examples all admit the Dirac-style formulations in terms of Hamiltonian 
vector fields. The integrability criterion and the first integrability algorithm formulated in Section 4 
can be adapted to this situation. 

Let Q be the configuration manifold of a mechanical system and let the dynamics of the system be 
represented by the union 

D=[j {im(X„|(7)} (460) 

aeA 

of a family of Hamiltonian vector fields 

X„:T*Q^TT*Q (461) 

generated by a family of Hamiltonians 

ffa:T*Q^R (462) 
and restricted to a submanifold C C T*(5 specified as 

C={p&l*Q; \jj,^A{p) = 0} , (463) 

where <^a are independent functions on T*(5 called primary constraints. The condition D C TC means 
that at points of C the vector fields X^. are tangent to C or that 

(d$A,^a)|C = (464) 

for each aeA and each function In view of 

the integrability criterion for the dynamics D assumes the form 

{F„,$a}|C = (466) 

for each aeA and each function 

Example 21. The dynamics of the system in Example 9 is integrable since it is the image of a 
Hamiltonian vector field. A 

Example 22. The dynamics in Example 10 is a Dirac system. It is the union of images of the family 
of Hamiltonian vector fields generated by the family 

Ha{q, q\p,Pj) = -^g'^ {Pi - eA,){pj - eAj) + e<f- a{p - e) (467) 

restricted to a constraint set C. The parameter a is an arbitrary function on T*Q. There is one 
primary constraint. It is the function 

^{q,q\p,Pj)=p-e. (468) 

The system is integrable since 

{H^,^}\C = 0. (469) 

A 
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Example 23. The dynamics in Example 11 is a generalized Dirac system. It is described by the 
family of Hamiltonian vector fields generated by the family 



Ha{q'',Px) = a - eA,){px - eA^) - mj (470) 

of Hamiltonians parametrized by a function a > on T*Q. There is one primary constraint 

Hq'',Px) = ^jg^^ipK - eA^){px - eAx) - m. (471) 
The integrability criterion is satisfied. 

▲ 

Example 24. The generalized Dirac system of Example 12 is described by Hamiltonian vector fields 
generated by the family 

H(a\a^){Q'',Px) (^^/g^^iPK - eA^){px - eAx) - mj +a^{p-e) (472) 

parametrized by functions > and on T*Q. There are two primary constraints; 



^i{q,q'',P,Px) = \I9''\Pk - eA^){px - eAx) - m (473) 

and 

*2(g,g^P,PA) =P-e. (474) 
Integrability criteria are again satisfied. ▲ 
Example 25. For the generalized Dirac system in Example 13 we use the family of Hamiltonians 

Ha{q'',Px) = ^9^^P.PX (475) 
depending on a function a > 0. There is one primary constraint 

*(9",Pa) =rt«PA. (476) 

The system is intcgrable. A 

Example 26. Let Q be the afiine space-time of special relativity with Cartesian coordinates (g**) 
and a constant Minkowski metric tensor (gnx)- We analyse the dynamics of two interacting relativistic 
particles [15] [16]. The configuration space is the product QxQ with coordinates (gf , (/2 )• Coordinates 
(9i,92>^f,92), {qi,q2,Pl,Pl), and {q^,q^,pl,pl,q^,qS,pl,pl) wiU be used in T(Q x Q), T*(Q x Q), 
and TT*{Q x Q) respectively. Masses of the particles are denoted by mi and m2. The interaction 
potential is a function F of a real positive argument. Relations 

g.xqiqi > 

g^xq2q2 > o 

g^V^x > (477) 

g'^VA > 

9Kx{q2 - 9f )(92 - 9i ) < 
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are assumed to be satisfied. Abbreviations 



1^2 - qi 



Pi 



(478) 



mi 



X2 - 



TO2 = \/m2^ + V{\\X2 - Xi\\) 

will be used. 

The dynamics of the particles is a generalized Dirac system described by Lagrange equations 



pi 



pI 



Pi 



pI 



W(l|g2-gi||) fm 
2||g2-gi|| V ^1 

Dmig2-giH) nvhW 

2||g2-?i|| V "^1 

mi 
m2 



mi 

m,2 

M 

m2 



gK\{q2 - ii) 
g^xiqi - 92) 



5«Ag2 



derived from the Lagrangian 

L{qi.q2Ai,q2)=mi \\qi\\ +m2||g2||. 
The Hamiltonian is the Morse family 

H{qx,q2,pl,pl,v^,v'^)=v^{\\pi\\-mi)+v'^{\\p2 
of functions of the variables > and > 0. Hamilton equations 

0^^(1192-9111) 



m2) 



Pi 



pI 



qi 



2ll'?2 - gill 

BV{\\q2-qi\\) 
2|k2-«i|| 



— + — 9.x{q2 - qt) 

mi 1712 J 

^ + 5«A(gi - 92) 
mi m2 / 



lb 



Til^ P\ 



«2 = ^^'^'f'i 

„i II 



2|| ™ 



nil 
m2 



are equivalent to the Lagrange equations (479). 

The dynamics is the union of the family of Hamiltonian vector fields 

^(ai,a^):T*(QxQ)^TT*(gxQ) 
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(479) 



(480) 



(481) 



(482) 



(483) 



generated by the family of Hamiltonians 

H(c.\c.^)iqi,q2,pl,pl) = a' {\\pi\\ - mi) + (Ib2|| - ma) (484) 
and restricted to a submanifold C CT*{Q x Q) described by two primary constraints: 

^i{Qi,Q2,Plpl) = \\Pi\\-mi (485) 

and 

^2{q^,q^,pl,pl) = \\p^-m2. (486) 
The parameters (a^,a^) are positive functions on T*{Q x Q). The Poisson brackets 

{i^(a^«^), = ^ %^i'i~''^!^M + plM - q^) (487) 



and 

{H,^^,.^,,^2]\C = - qr) 

zmi 7712 1192 — qi\\ 

indicate that the dynamics is not integrable unless the potential V is constant. 



(488) 



The generalized Dirac system in Example 26 is the only case of a non integrable dynamics known to 
us. We will extract the integrable part of this system by applying different versions of the extraction 
algorithm. 

Example 27. We apply the first algorithm of Section 4 to Hamilton equations (482) and constraint 
set C described by 



IIp^II = "T-i 

11/11 =m2 



(489) 



We obtain equations 



Ib^ll = "^1 
Ib^ll = m2 

9^VA = "^.["'^^"^if g.A(g2- - q^M - q^) (490) 

2||g2 - gill 

„kA„2 .2 _ D^(||g2-gl||) . « A •A^ 

9 PkP\ - — i|— fi'KA(?2 - ^1 j(«2 - «1 j 



for TC and equations 



1 _ Dl^(||g2-gi|i; <- , X \^ 

Pk = —^i n — — + — aK\{q2 - qi ) 

2||<?2-gi|| V"''i ™2/ 

.2 T)V{\\q2-qi\\) f ,x y^. 

pI = on „ II [^ + ^] 9nxiq^ - 92) 



V 



A„l 



' A 2 (491) 



^2 = J^^9 p. 

\\p'\\=rni 
\\p'\\=rn2 

o = {pi+pl)iq2-q1) 

v^>0 
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for the intersection = D (1 TC if the potential V is not constant. For the new constraint set 
C = tq{D ) we have equations 

lb^ll=mi 

lb' II =^2 (492) 
= {pl+pl){q^-qt) 



and equations 



\\p'\\=rn, 
lb'll=^2 

g'^Vji = ^^^9.M - qDi<i^ - q^) (493) 
9^VA = ''\!l^''-'f 9.M-q?M~q^) 

Am - qi\\ 

= m+plM - qi) + (pI+pIM - q'i) 



for Tc'. The equations 



pI 


D^(lk2-9i||) 

21192-9111 


\mi 


„2 N 
m2j 


) 9K\{q2 




pI 


_ Dl/(||g2-gi||) 
2|k2-9i|| 


fv^ 
Vmi 


y-2 \ 


] 9^x{qi 


-q^) 


ql 


= iiP'ii'""^ 












v'^ k\ 2 










lb' II 


= mi 










Ib'll 


= fn2 











(494) 

= m2 

pI + pI^^ 

0^{pl+plM-q^i) 

^9''\pI+pI)p\ = ^9^\pI+pI)pI 
nil ^2 



2 1 2 

for = D n TC of D are the last step in the algorithm. We have excluded from D the case 
p\+p\ = 0. The system D = D will be shown to be integrable. A 



Example 28. We apply the second algorithm of Section 4 to the Lagrange equations (479). The 
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prolongation of the Lagrange equations is the system of second order equations 

, _ DV{\\q2-qi\\) f\M , 119211 \ ^ .A X, 



2||(?2 - gill V "^1 '^2 ' 

. BV{\\q2-qi\\) fUiW , ||g2|| ^ „ ,x A^ 



21192 -qi\\ \mi 1712 J 

1 ^1 -A 
p. = ^.«A.l 

2 ^2 .A 

^'^ = M^^"^' (495) 

pl = f\qi,q2,<i'i, 41^1,92) 
pI = f{ql,qlqlqliilii2) 

^« = - miiir-'^'ii ^-^^-' - ^^^^^^^ - w + ^^^^^ ('^' - 

.2 W(l|g2-glH) ^w., 92 ^_ Z^rA 



qW^] 


1 


\m) 




q^q'2] 


1 € 


\m) 


' 119211 



The exact form of the functions /^(^f , 92 ) 9]*) 92 ) 9i ' 92 ) ^^'^ /^(9i ; 92 ) 9f ) 92 ) 9i ' 92 ) is of no interest 
since equations for p]. and can not impose restrictions on first derivatives. Projection in T(Q x Q) 
eUminates these equations and also the components of the last two equations orthogonal to {qf) and 
(92 ) respectively. The resulting first order equations 

, _ Dl-(||92-9l||) ni9l|| , ||92||\ „ .A .A^ 



P. = o I ^ + ^ 5«a(9^ - qt) 

2||92 - 91 II V "ii "^2 / 

, Dnil92-9l||) /II91II , l|92|iy^ .A „A^ 

- _2||92-9i|| I W + ^'^^^^^ " 

1 JTli .A 

Il9i|| 

2 ^2 -A 

P. = ^5.A92 

2||92-9i|| mi 

.2.«_ Dmi92-9l||) ||92||^ .A M.. 

^'^'^ " " 2||92-9i|| ~ "'^'^ 



(496) 



are equivalent to the simplified equations 



pI = 


Dmi92-9i||) 


2||92 - 91 II 


pI = 


Dmi92-9i||) 


2||92-9i|| 


Pi= 


m,i .A 

Ikiill 


pI- 


m2 .A 

II • II 5kA92 

119211 


= 


{pI+pIM- 



''^^"+^).«a(92^-9^) 



mi m2 

P + P)..a(9N92^) 
mi m2 / 



(497) 



These equations are not yet integrable. They were falsely declared the integrable part of the dynamics 
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in [8]. The prolongation of equations (497) projected in TT (Q x Q) results in equations 

+ 9Kx{q2 - Qi) 





Dy(||g2-9i||) 


2||g2 - 91 II 


Pk — 


DV(|k2-<Zi||) 


2|k2-gi|| 


Pi = 


mi 
Il9i|| 


pI = 


TO2 .A 
II • II 9ii\l2 

Ik2|| 


= 


{pI+pI){i2- 


= 


ipl+pl){Q2- 



mi m2 
mi m2 



+ SkaIO'i - 92. 



(498) 



Prom 

{Pi+Plm - Qi) = ^^9^\pi +pI)pI - ^%^9^\pi+pl)pi (499) 
m2 mi 

we see that the last equation in (498) imposes a synchronization relation between parametrizations of 
the world lines of the particles unless P^+pf^ — 0. The case p}^ + pf^ = seems to be too restrictive 
to be of interest. We will exclude this case. The integrability algorithm applied to equations (498) 
produces no further restrictions. A 

Example 29. Although the dynamics of two relativistic particles is not a Dirac system it admits the 
Dirac-stylc representation in terms of constrained Hamiltonian vector fields. This representation can be 
used to simplify the algorithm in Example 27 and to prove the integrability of the resulting equations. 
The simplified algorithm is an adaptation of Dirac's original algorithm [1]. Since the Poisson brackets 
(487) and (488) do not vanish on C we restrict the set C by adding the secondary constraint 

*(9f , q^,pl,pl) = {pi +pIM - Qi)- (500) 

The resulting system is now the union of images of Hamiltonian vector fields generated by the 
family (484) restricted to the constraint set C CT*{Q x Q) satisfying the equations 

$1 = 0, $2 = 0, * = 0. (501) 

The system is not yet integrable. In the second step we require the vanishing of the Poisson bracket 



,,1 ^,2 



on the new constraint set. We will exclude the case + = and impose the condition 



^9^\pI+pI)p\ - ^9^\pI+pI)pI = (503) 
mi m2 



on the functions (a^ , a^). The resulting system D is the union of images of constrained Hamiltonian 
vector fields generated by the family (484) with the positive functions (Q:^,a^) satisfying the condition 
(503). The vector fields are restricted to the constraint set C . The system is exactly the system 
described by equations (494) and (498). The system is integrable since 

{F(„i,„2),$i}|c'=0, (504) 

= 0, (505) 

and 

{F(„i,„2),*}|c'= 0. (506) 

▲ 
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